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COMVF.RSION  FACTORS,  U.S.  CUSTONF\RY  TO  MI'.TRIC  (SIJ 
UNITS  OF  MHASURFiMFiNT 


U.S.  customary  units  of  measurement  used  in  tliis  report  can  he  converted 
to  metric  (SI)  units  as  follows; 


Multiply  by  To  obtain 


inches 

25.4 

mi  1 1 imeters 

2.54 

centimeters 

square  inches 

b.452 

square  centimeters 

cubic  inches 

lb.  39 

cubic  centimeters 

feet 

30.  39 

centimeters 

0. 3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

ki lometers 

square  miles 

259.0 

hectares 

knots 

1.8532 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  X 10 

kilograms  per  square  centimeter 

ounces 

28.  35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0. 1745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins* 

^To  obtain  Celsius 

(C)  temperature 

readings  from  Farcnheit  (F)  readings, 

use  formula:  C = 

(5/9)  (F-32). 

To  obtain  Kelvin 

(K)  readings,  ase 

formula:  K = (5/9)  (F  -32)  273.  15 
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turbulent  velocity  fluctuations 

k wave  numlier 

L amplitude  of  flume  oscillation 

1 length  of  surface  wave 

Ig  length  scale  for  oscillating  and  unidirectional  flow 

M exponential  decay  rate  of  sediment  concentration  of  oscillating 

flow 
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R hydraulic  radius 

r horizontal  oscillating  velocity 
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!^noma8  C.  MaaDonald 

I.  INTRODUCTION 

SeJiinent  transport  by  waves  approaching  the  shore  has  been  analyzed 
in  two  ways.  Inshore  of  the  breaker  zone,  the  extremely  complex  flow 
patterns  and  turbulence  resulting  from  the  breaking  waves  have  necessi- 
tated only  a qualitative  approach  to  sediment  transport  with  quantita- 
tive estimates  based  on  field  measurements.  Offshore  of  the  breaker 
zone  in  relatively  deep  water,  the  problem  of  sediment  transport  can  be 
approached  in  a more  theoretical  manner.  In  this  zone,  sediment  trans- 
port studies  are  simplified  because  there  is  no  turbulence  in  the  flow- 
field  from  breaking  waves  and  the  flow  condition  near  the  ocean  bottom 
can  be  better  estimated  from  linear  wave  theory.  This  report  concerns 
sediment  transport  offshore  of  the  breaker  zone. 

Laboratory  and  field  observations  indicate  that,  as  in  unidirection- 
al flow,  sediment  transport  at  the  ocean  bottom  offshore  of  the  breaker 
zone  is  of  two  different  types--bedload  and  suspended  load.  The  dis- 
tinguishing feature  between  these  two  tyjjes  of  transport  is  that  in 
suspended  transport  the  entire  weight  of  the  sediment  is  continuously 
supported  by  the  fluid;  whereas,  in  bedload  transport  the  sediment  rolls, 
skips,  and  jumps  along  the  bed  and  therefore  its  weight  is  partially 
supported  by  the  stationary  bed.  For  moving  sediment  to  be  supported 
by  the  bed  means  that  the  regime  of  bedload  transport  is  contained  in  a 
thin  layer  adjacent  to  the  stationary  bed,  two-grain  diameters  thick  as 
proposed  by  Einstein  (1950).  The  majority  of  the  sediment  in  motion  in 
this  area  is  bedload  and  thus  most  research  has  concentrated  on  the 
turbulent  boundary  layer  and  the  oscillatory  hedload  rate  due  to  wave 
action  (Li,  1954;  Manohar,  1955;  Kalkanis,  1957,  1964;  Abou-Seida, 

1965).  Sufficient  advances  in  the  theory  of  bedload  movement  in  an 
oscillating  flow  have  warranted  studying  the  suspended  load  to  deter- 
mine: (a)  Approximately,  what  percentage  of  offshore  movement  is  duo 

to  suspended  load,  i.e.,  a second-order  approximation  to  total  trans- 
port; and  (b)  if  any  of  the  now-predicted  bedload  is  partially  suspended 
load. 

This  investigation  develops,  from  an  empirical  approach,  a method 
for  predicting  the  distribution  of  suspended-sediment  concentration 
based  on  the  hydraulic  flow  conditions;  i.e.,  surface  wave  amplitude 
and  period,  depth,  sediment  characteristics,  and  bottom  roughness  con- 
ditions. Although  only  one  bottom  roughness  was  studied,  the  resulting 
method  is  general  enough  to  be  extended  to  other  rougfiness  conditions 
l)y  additional  experimentation.  Ilic  suspended  distributions,  when  ased 
in  conjunction  with  the  bedload  function  of  Kalkanis  (1964),  should 
give  a better  approximation  of  the  total  sediment  transport. 


As  in  Kalkaiiis'  bedload  function,  the  approadi  to  the  suspended 
load  is  based  on  many  of  the  same  principles  proposed  l)y  Hinstein  flDSO) 
in  his  theory  of  bedload  and  suspended- load  transport  in  unidirectional 
flow.  Analysis  of  suspended  load  in  oscillating  flow  is  more  complicated 
than  that  of  unidirectional  flow  because  of  two  factors.  First,  in  uni- 
directional open  channel  flow  the  entire  depth  of  flow  is  turbulent  and 
the  relatively  liigh  turbulent  velocity  fluctuations  allow  the  sediment 
exchange  coefficient  to  be  approximated  by  the  nnimentum  exchange  coeffi- 
cient which  can  be  obtained  from  the  shear-stress  distribution.  In 
oscillatory  flow  this  is  not  possible.  Both  laboratory  and  field  obser- 
vations indicate  that  suspension  of  sediment  occurs  to  depths  considerably 
abov'e  the  boundary  layer  in  an  area  where  the  shear  stresses  duo  to  the 
oscillating  motion  are  extremely  small  and  difficult  to  measure.  Although 
it  may  be  possible  to  express  a sediment  exchange  coefficient  in  the 
boundary  layer  as  a function  of  the  moan  sliear  stress,  this  would  not 
provide  a means  of  estimating  the  sediment  exchange  coefficient  above 
the  boundary  layer.  Tlierofore,  a sediment  exchange  coefficient  which  is 
not  based  on  a shear-stress  distribution  must  be  found. 


The  second  factor  concerns  the  magnitude  of  the  turbulent  velocity 
fluctuations.  Offshore  of  the  breaker  zone  where  the  flow  velocities 
near  the  bed  are  low,  only  a small  part  of  the  wave  energy  is  dissipated 
by  friction  at  the  boundary.  The  remaining  wave  energy  is  lost  inshore 
by  the  breaking  waves.  Because  of  the  relatively  low  intensity  of  tur- 
bulence in  tliis  offshore  area,  the  vertical  velocity  fluctuations  are  of 
the  same  order  of  magnitude  as  the  settling  velocity  of  the  sediment. 

Under  these  conditions,  the  sediment  excliange  coefficient  is  highly 
dependent  on  the  sediment-settling  velocity.  Therefore,  measured  dis- 
tributions of  turbulent  velocity  fluctuations  and  sediment  concentrations 
mast  be  used  in  analyzing  the  upward  turbulent  flux  and  downward  turbulent 
and  gravitational  flux  for  each  sediment-settling  velocity. 

To  obtain  a relationship  between  sediment  suspension  and  flow  hydrau- 
lics in  an  oscillating  flow,  concentration  distributions  for  various  flow- 
conditions  must  be  measured.  The  sediment  exchange  coefficient  is  deter- 
mined from  these  measurements.  Next,  the  turbulent  velocity  fluctuation 
distribution  with  time  at  a constant  elevation  and  its  distribution  with 
elevation  mast  be  measured.  This  measurement  will  yield  the  information 
necessary  to  describe  the  fluid  exchange.  The  distribution  with  eleva- 
tion will  yield  a velocity  scale,  one  of  the  two  variables  composing  the 
sediment  exchange  coefficient.  From  the  sediment  exchange  coefficient 
and  the  velocity  scale,  the  second  variable  (the  length  scale  or  its 
associated  time  scale)  can  be  calculated.  Knowledge  of  these  fundamental 
variables  of  suspension  as  a function  of  the  flow  hydraulics  stiould  lead 
to  a practical  method  of  estimating  the  suspended- load  dist ri!>ut i on , and 
indicate  the  important  variables  of  suspension  in  an  oscillating  flow. 

ihe  only  otlier  retiuirement  for  a solution  to  tlie  suspended  load  is  a 
knowledge  of  a base  concentration  as  a function  of  flow  hvdraulics.  Ihe 
base  concentration  is  determined  from  Kalkanis'  (IPOd)  beiiload  theory, 
lieing  the  concentration  at  tlie  top  of  tlie  bedload  layer. 
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Combination  of  the  t)edload  and  suspended  load  will  predict  the  total 
amount  of  sediment  in  motion  under  specified  wave  and  l)oundar>'  conditions. 
By  superimposing  a constant  unidirectional  flow,  such  as  tlie  mass  trans- 
port or  a coastal  current,  the  amount  of  sediment  transport  can  be 
est imated. 


In  this  investigation,  as  in  many  research  projects,  the  experimental 
researcli  preceded  the  development  of  a suitaiile  method  of  describing  sedi- 
ment suspension.  For  this  reason  a description  of  tiie  experiments  and 
their  results  will  be  given  first,  followed  by  a discussion  of  how  the 
results  can  be  used  to  predict  the  sedimoit  suspension  load. 

11.  CO\Ci:.NTR..\TION  DISTRIBUTIONS 
1 . Hxperimental  Apparatas. 

E.xperiments  and  observations  indicate  tliat  near  the  ocean  bottom  off- 
shore of  the  breaker  zone  in  relatively  deep  water,  sediment  is  held  in 
suspension.  Tliis  is  due  to  the  turbulence  resulting  from  the  dissipation 
of  wave  energy  on  the  rough  ocean  bed.  For  waves  with  a small  surface 
slope,  where  3/3x  « 3/3y,  the  fluid  motion  c;iri  he  .ipproxi  m.Jted  from 
linear  wave  theory.  Tlie  equations  describing  the  horizontal  and  vertical 
displacement  of  a fluid  particle  are  given  by  the  expressions  from  iamb 
fl9  32)  : 


X = 


{cosli  [k  (y*d  ) J/s  i nh  (k  d)}  cos  (k  x - w t) 


(1) 


A = {sinh  [k (y  + d) ]/sinh (k  d|}  sinfk  x - w t]  (2) 

where 

d = the  water  depth 

y = the  distance  from  the  mean  water  surface  me.isured  negatively 
downward 

11  = the  crest-to-trough  wave  height 
k = 2tt/1 

1 = the  length  of  the  surface  wave 
w = 2it/T 

T = the  period  of  the  surface  wave 
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nu'  coirespojuling  velocity  components  are  obtained  by  differentiation 
with  respect  to  time  of  tiie  above  equation^  to  give: 

u = w {cosh  [k  (y+d)  J/s  inh  fk  dj}  sin(k  x - w t)  '’3) 

r = ("4Hj  w (sinli  [k  (y+d)  ]/s  inli  (k  d)  } cos  (k  x - w tj  (4) 

i'rom  these  equations  it  is  evident  that  the  vertical  component  of  the 
displacement  and  velocity  liecomes  smaller  as  the  distance  from  the  sur- 
face increases.  At  the  bottom,  whore  y = -d,  the  motion  degenerates  into 
a simple  harmonic  oscillation  in  the  x-direction.  It  is  this  horizontal 
harmonic  oscillation  which  is  of  first-order  imjiortance  in  producing  tur- 
bulence and  suspension  of  sediment.  Tlie  equations  also  indicate  that  the 
magnitudes  of  the  horizontal  displacement  and  velocity  change  slightly 
with  depth  near  the  ocean  bottom  and  can  therefore  be  considered  constant 
in  tile  region  of  sediment  saspension.  For  example,  consider  a S- foot-high 
wave  with  a wavelengtli  of  l.SO  feet  and  a period  of  10  seconds  in  a water 
depth  of  50  feet,  liquation  (3)  indicates  the  horizontal  velocity  at 
y = -d  is  11  = 0.393  sin(kx-wt)  feet  per  second  and  that  for  y = -d  + 5 
(5  feet  above  tlie  bed),  the  horizontal  velocity  is  u = 0.401  sin(kx-wt). 
ilie  sediment  suspension  measurements  which  will  be  discussed  later  indi- 
cate that  the  regime  of  measurable  sediment  concentrations  is  well  with- 
in the  bottom  5 feet  of  water  depth  for  this  typical  wave  condition. 
(Measurable  sediment  concentrations  were  actually  found  within  1 foot 
of 'the  bed.)  Ilie  change  in  horizontal  velocity  in  the  bottom  5 feet  of 
water  depth  is  only  2 percent.  Tlierefore,  to  design  an  experimental 
apparatus  to  simulate  the  flow  conditions  near  the  ocean  floor,  the  hori- 
zontal flow  velocity  above  the  boundary  layer  can  be  considered  constant 
and  equal  to  the  value  given  by  equation  (3)  for  y = -d.  Equation  (4) 
indicates  that  for  the  wave  condition  discussed  above  the  vertical  flow 
velocity  at  y = -d  is  zero,  and  at  an  elevation  of  5 feet  above  the  beil 
is  0.082  cos  (kx-wt) . Tlie  fact  that  the  vertical  velocity  only  increases 
slightly  in  the  5 feet  above  the  bod  and  that  its  motion  is  symmetrical 
suggests  that  the  vertical  velocity  has  no  effect  on  the  suspension  of 
sediment. 

With  the  above  assumptions  of  a constant  horizontal  oscillating 
velocity  and  zero  vertical  oscillating  velocity  above  the  boundar>’  layer, 
the  turbulent  flow  conditions  which  exist  near  the  ocean  floor  arc  easily 
approximated.  By  superimposing  a constant  velocity  equal  to  that  given 
by  equation  (3)  for  y = -d  but  opposite  in  sign,  there  would  lie  no  motion 
in  the  fluid  above  the  boundarx'  layer,  the  distribution  of  velocities  in 
the  boundary  layer  would  l)c  inverted,  and  the  bed  would  be  oscillating  at 
the  simple  harmonic  given  by  e(|uation  (3)  for  y = -d. 

It  is  now  only  necessary  that  the  physical  apparatus  wh  i cli  du]ilicates 
these  flow  conditions  contains  a water  depth  greater  than  tlie  thickness 
of  the  boundary  layer.  Kalkanis  (1957)  made  velocity  distribution  measure- 
ments in  a flume  containing  a still  body  of  water  with  an  oscillating 


rough  bod.  Kalkanis  (19ti4)  showed  tliat  the  velocity  distribution  in  an 
oscillating  flow  with  the  same  assumptions  as  described  above  can  be 
approximated  by: 

u/u,T  = [1  + fj(Y)  - 2f^(Y)  cos  f2(Y)]  - sin(wt  + 6)  (5) 

where 

fl(Y)  = exp(Y- loVaSD)  , (6) 

fzCY)  = 0.5(BY)0-67  , (7) 

and  sin(wt  + 0)  describes  the  variance  of  the  velocity  witli  time  and 
phase  angle.  In  equations  (6)  and  (7),  Y is  the  elevation  above  the 
bed  measured  jiositively  upwards;  a is  the  amplitude  of  bed  oscillation; 

6 = (w/2v)  ■^;  V is  the  kinematic  viscosity;  and  D is  the  representa- 
tive roughness  diameter.  For  the  flume  roughness  conditions  used  in  this 
investigation  (D  = 0.05  foot),  and  the  flow  velocities  studied  [0.67  foot 
< a < 2.0  feet  and  2.0  seconds  < T < 15.0  seconds) , equation  (5)  indicates 
the  boundary  layer  thickness  is  no  more  than  a few  millimeters  thick. 

To  approximate  the  shear  stresses  and  therefore  the  turbulence  condi- 
tions near  the  ocean  bed  for  a given  water  depth,  wave  period,  and  an^ili- 
tude,  a flume  having  a moving  bed  under  a still  body  of  water  was  used. 

The  frame  of  reference,  which  describes  the  prototype  fluid  motion  under 
the  above  assumed  conditions,  is  then  moving  with  the  bed.  It  is  only 
necessar>'  to  oscillate  the  bed  in  the  harmonic  motion  described  by  the 
linear  wave  theory  for  y = -d.  The  swing  flume  used  in  this  investiga- 
tion duplicates,  on  a one-to-one  scale,  these  conditions. 

The  swing  flume  is  shown  in  Figure  1(a).  The  flume  bed  is  shaped  to 
an  arc  segment  of  a circle  with  an  8.92-foot  radius,  a 13.. 5.5- foot  chord 
length,  and  a 12-inch  width.  The  flume,  suspended  from  the  ceiling  of 
tlie  laboratorv’,  is  free  to  rotate  about  its  center  of  curvature.  Tlie 
flume  is  oscillated  about  its  center  position  by  a 1.5-horsepower 
variable-speed  motor  connected  to  a drive  wheel  with  an  eccentric  arm. 

The  eccentric  arm  is  connected  by  a 10-foot  connecting  rod  to  a linkage 
fixed  to  the  flume  bottom.  The  linkage  at  the  flume  is  adjustable  to 
allow  correction  of  the  asymmetry  of  motion  which  would  result  from  a 
change  in  eccentricity.  Variations  in  eccentricity  and  motor  speed  allow 
the  amplitude  and  freuuency  of  oscillation  to  be  varied  over  a wide  range 
of  prototype  wave  conditions. 

Within  the  flume  is  a stationary  horizontal  board,  slightly  less  than 
12  inches  wide,  8 feet  long,  and  at  an  elevation  of  12  indies  above  the 
lowest  point  of  the  curved  liottom.  Tlie  board,  which  is  separately  su]i- 
ported  from  the  ceiling  of  the  building  and  is  not  connected  to  the  flume, 
suppresses  any  st.inding  surface  waves  in  the  flume  caused  by  the  flume 
motion.  riu'  fluid  at  this  elevation  must  remain  stationarv'  to  conform 
to  the  flow  comlitions  described  above. 


I 


4 


15 


^ UNa-ASSlIlID 

SECUOITV  CLASSIFICATION  OF  THIS  PAGEITFhF'i  Oaif  Fnl0ima> 


Asymmetric  End  Roughness- 

Two-Dimensionol  Roughness  Elements 
/ rPoint  of  Suspension 


U-  A 

V '-e"  ot; 


b.  Cross  Section 
of  Roughness 
Elements 


0.9"  0.9" 

i L 


c.  Cross  Section  of 
Asymmetric  End 
Roughness 


^ l3'-4"- 

0 Elevotion  View  of  Entire  Flume 


Gloss  Window 


Horizontal  Board  to 
Suppress  Standing 
Surface  Waves 


Figure  1.  Swing  flume. 
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The  shape  of  the  artificial  bed  rougiiness  (Fig.  l,(i)  was  deternined 
by  experiments.  A large  quantity  of  the  artificial  jilastic  sediment  used 
in  concentration  measurenents  was  put  into  tlie  smootli-bot tomed  flume. 

Tile  flume  was  oscillated  at  various  amplitudes  and  periods  covering  the 
range  of  flow  conditions  to  he  studied.  After  the  flume  was  oscillated 
at  a constant  rate  for  a period  long  eno'ngli  to  establish  a natural  bed 
shape,  the  flume  was  stoiiped  and  the  bed  dunes  were  ne.asurc'd.  iiie  bed 
shape  was  found  to  be  approximately  sinusoidal  in  the  cross  section  under 
all  flow  conditions.  The  mean  wavelength  of  the  bed  shape  was  .A. 5 inches 
with  a range  of  4..S  to  8.0  inches;  the  mean  wavel  ength- to-depth  ratio  was 
8.0  with  a range  of  to  12.  The  fixed  artifici.al  roughness  used  apjirox- 
imates  this  shape.  Tlie  artificial  dunes  were  constructed  of  wood  and 
fastened  to  a flexible  sheet  of  plastic.  Natural  sediment  with  a mean 
diiuneter  of  0.3  millimeter  was  glued  to  both  the  dunes  and  [ilastic.  ihe 
plastic  sheet  was  fixed  to  the  flume  bottom,  covering  the  central  (>.53 
feet  of  the  arc. 

Tlie  asvTinetric  end  roughnesses  shown  in  Figure  1(c)  were  used  to 
eliminate  secondary  currents  in  the  central  measuring  section  of  the 
flume.  Proper  placement  of  tlie  asymmetric  roughness  elements  depends  on 
the  flow  conditions  of  the  flume,  ihe  optimum  placement  of  the  rougiiness 
elements  for  a given  flow  condition  was  determined  by  dropping  iiotassium 
dichronate  crystals  through  holes  in  the  horizontal  wave  suppressent  board 
and  observing  the  movement  of  the  dye  streaks.  IvTien  no  transi'orsc  move- 
ment of  the  dye  stre.aks  in  the  central  part  of  the  flume  were  observed, 
the  roughness  elements  were  considered  to  be  in  the  Ojitimum  location. 

A lightweight,  black  plastic  material  with  specific  gravity  of  1.2.S, 
which  had  been  crushed  and  sieved,  was  used  as  artificial  sediment  in  the 
experiments.  The  grain  diameter  of  the  sediments  was  uniform,  bracketed 
by  two  consecutive  sieve  sizes  of  approximately  the  same  diameter  as  the 
natural  sediment  glued  to  the  flume  bed.  Only  a small  quantity  of  the 
sediment  was  used  in  tlio  flume  during  an  experiment  in  order  to  limit 
the  deposition  of  sediment  wliich  would  alter  the  flume  bottom  geometry. 

Wlten  tlie  sediment  was  first  put  into  the  flume  it  was  found  that  air 
bubbles  adhered  to  the  sediment  particles,  thereby  dianging  its  settling 
velocity.  To  eliminate  tliis  buoyancy  effect,  deaerated  water  was  used. 

Ihe  water  was  deaerated  in  a 60-cub ic  foot- capacity  tank  located  on  tlie 
wall  of  the  laboratory  at  an  elevation  above  the  swing  flume,  heated  by 
a 5-kilowatt  immersion  heater  to  a temperature  of  90*^  Fahrenheit,  and 
then  cooled  to  room  temperature.  Tlie  deaerated  water  was  transported  to 
the  swing  flume  by  gravity  through  a hose  to  reduce  air  entrainment. 

Tlie  optical  concentration  meter  used  in  the  experiments  was  dovoloiicd 
by  Das  (1971)  for  measuring  in  pitu  concentrations  in  laboratory  flumes. 

Tlie  equipment  consists  of  a light  source,  a beam  collimator,  a receiving 
unit,  and  necessary  recording  units.  A collimated  beam  of  light,  8..5- 
raillimeter  average  diameter,  is  projected  through  the  glass  walls  of  the 
flume  to  a duo-photodiode  mounted  on  the  opposite  siile  of  the  flume. 
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I'hc  ptiotOL' U'lt  ri  c cell  produces  ;i  signal  wliich  is  proportional  to  the 
amount  of  lig/it  received  which,  in  turn,  is  proportional  to  the  amount 
of  light  blocked  out  by  suspended  sediment.  Tlie  signal  from  the  photo- 
electric cell  is  transmitted,  amplified,  and  recorded  on  an  analog  paper 
chart  recorder  and  an  analog-to-di gital  data  acquisition  system  set  to 
sample  at  the  rate  of  58  samples  per  second.  'Ilie  liglit  source  and 
receiver  are  shown  in  l-igure  2. 

Hie  light  source  and  receiver  are  inoimted  on  a rigid,  aluminum  yoke 
support  (fig.  held  by  friction  brackets;  by  loosening  four  thumliscrews 
the  elevation  of  the  sup[)ort  can  be  varieil  from  below  the  flume  bottom  to 
above  the  flume  top  wiiile  maintaining  precise  alinement  of  the  optical 
equipment.  Tli  i s fle.xibility  allows  calibration  of  the  optical  equipment 
and  flume  concentration  measurements  to  bo  made  without  removing  the 
equitnuent  from  the  suf>ports.  Itie  )oke  support  and  brackets  are  pinned 
to  the  Laboratory  ceiling  on  the  same  a.xis  as  the  flume  so  that  the 
optical  equipraint  can  swing  with  the  flume  or  be  held  stationary  in 
space . 

“•  h-xperimental  I'rocedure. 

a.  Settling  Velocity  of  the  Sediment.  Hie  sediment  used  in  the  first 

series  of  experiments  was  the  material  whicli  passed  the  0 . 4yS-mi  1 1 i meter 
sieve  and  was  retained  on  the  0.417-millimeter  sieve.  The  fall  velocity 
of  this  sediment  was  me.asured  in  a 2- inch-d i ameter  glass  cylinder  filled 
with  deaerated  water  at  a temperature  of  72°  Falirenheit.  llie  time  rc- 
ciuired  for  220  particles  chosen  at  random  to  fall  8.50  inches  was  measured 
with  a stopwatch.  Tlie  average  settling  velocity  the  range  of 

velocities,  and  the  standard  deviation  were  then  calculated  to  be  0.055, 
0.0213  to  0.0532,  and  0.0059  foot  per  second,  respectively. 

b.  Calibration  of  Optical  Concentration  Meter.  Calibration  of  the 
optical  equipment  was  done  in  a 0.5-  by  0.5-  by  1.0-foot  clear,  plastic 
calibration  tank  placed  on  the  top  of  the  swing  flume.  ilie  1.0-foot 
dimension  of  the  tank  was  positioned  parallel  to  tlie  axis  of  tlie  light 
source  and  receiver;  i.e.,  the  same  width  as  the  swing  flume.  Tlie  tank 
was  then  filled  with  a measured  quantity  of  deaerated  water.  A small, 
measured  amoimt  of  cleaned  sediment  was  added  to  the  tank  and  the 
sediment-water  mixture  was  stirred  mechanically  to  give  a luiiform  sus- 
pended concentration.  Tlie  collimated  light  beam  and  receiver  were  posi- 
tioned and  a 5-secoiul  record  of  the  voltage  from  the  receiver  was  recorded 
on  tlie  analog  chart  recorder  and  magnetic  tape.  Uniformity  of  concentra- 
tion in  the  tank  was  checked  by  measuring  the  fraction  of  light  blocked 

by  the  sediment  at  various  locations  in  the  cross  section  of  the  tank. 

Ihe  mechanical  stirrer  was  stopped  and  the  sediment  allowed  to  settle. 
Records  were  tlien  made  of  tlie  voltage  from  the  liglit  beam  passing  through 
clear  water  and  the  voltage  of  the  amliient  light.  Tlie  fraction  of  light 
passed  was  then  calculated  as  the  ratio  of  the  voltage  with  sediment  in 
suspension  to  the  difference  between  the  voltages  of  the  beam  through 
clear  water  and  the  ambient  light.  Measured  amounts  of  scdimt'iit  wore 
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adJod  to  tlic  calibration  tank  and  the  procedure  repeated  until  the  range 
of  tlie  optical  ei|uipinent  was  covered.  Measurable  concentrations  ranged 
from  0.1  to  2..")  grams  per  liter,  'llie  logarithm  of  the  fraction  of  ligiit 
passed  versus  tlie  concentrat ion  in  grams  per  liter  was  plotted  to  give 
the  calibration  curve  shown  in  figure  -1.  The  least  scpiares,  best  fit 
eciuation  for  this  curve  is: 

C = (-O.-SSa)  lng[T/(T^  - T')\  , (8J 


where 

C = concentration  (in  grams  jier  liter  I 

I = mean  voltage  from  the  receiver  when  the  light  beam  is  passed 
through  the  sediment-water  mi.xture 

= mean  voltage  from  tiie  receiver  when  the  light  beam  is  passed 
through  clear  water 

1'  = mean  voltage  from  tlie  receiver  when  tiie  light  beam  is  turned  off 

The  calibration  curve,  which  was  checked  periodically  during  the  experi- 
ments, did  not  change. 

In  setting  up  the  calibration  experiments,  it  was  found  that  a change 
in  the  focus  of  the  optical  eiiuipmcnt  would  cliange  the  resulting  calibra- 
tion curve.  To  eliminate  this  problem  a brace  was  made  to  hold  the  equip- 
ment in  focus.  A check  of  the  focus  was  made  during  each  concentration 
measurement  by  a wire-screen  filter  which,  when  placed  in  the  liglit  beam, 
blocked  out  a known  and  constant  amount  of  light.  If  a change  in  focus 
was  detected  by  the  filter  measurements,  the  concentration  measurements 
were  not  used. 

c.  Concentration  Distrihution  Measurements  For  concentration  dis- 
tribution measurements  in  the  flume,  a datum  elevation  of  the  optical 
equipment  had  to  he  established.  The  mean  elevation  of  the  crest  of  the 
artificial  dunes  wa.s  used  as  the  base  elevation  and  was  determined  by  a 
scale  fixed  to  the  flume  and  a pointer  fixed  to  the  optical  eipiipment 
brace. 

Tlie  desired  period  and  eccentricity  were  selected  and  the  flume  link- 
age adjusted  to  give  a symmetric  oscillation.  The  flume  was  filled  with 
deaerated  water  to  the  elevation  of  the  wave  suppressent  board  (12  inclics 
above  the  lowest  point  of  the  bed),  and  the  asymmetric  end  roughnesses 
were  positioned.  Depending  on  the  flow  conditions  chosen,  100  to  .'iOO 
grams  of  sediment  was  cleaned  and  deposited  in  the  flume.  Ibe  flume  was 
then  oscillated  until  the  distribution  of  sediment  in  the  flume  was  at 
equilibrium.  Tlie  equilibrium  condition  was  determined  by  periodically 
measuring  tlie  concentration  of  sediment  at  a fixed  point  in  space  until 
the  concentration  did  not  change  witli  time. 
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'llie  ooiK-ent  rat  ion  J Lst  r i luiT  i on  in  the  vertical  vvas  tlien  measured  for 
the  flow  condition  selected.  Again,  for  cacti  elevation,  voltage  records 
were  made  for  the  sediment- laden  water,  clear  water,  ambient  liglit,  and 
clear  water  plus  the  filter.  If  tlie  optical  equipment  was  held  stationary 
in  sjiace,  the  records  for  clear  water  and  clear  water  plus  filter  were 
made  while  the  flume  was  oscillated  at  sucli  a long  period  tliat  no  sediniert 
was  in  suspension.  llie  time  interval  of  the  record  for  sediment- laden 
water,  clear  water,  and  clear  water  plus  filter  was  always  an  integer 
multiple  of  the  flume  oscillation  period.  Tltis  procedure  automatically 
allowed  any  irregularities  in  the  transmissibi lity  of  the  flume  windows 
to  be  compensated  for  when  calculating  the  fraction  of  liglit  passed. 

'Hie  output  signal  of  the  photoelectric  cell,  which  was  recorded  on 
magnetic  tape  at  the  rate  of  58  samples  per  second,  was  not  constant 
with  time,  due  to  instantaneous  concentration  fluctuations.  Because  the 
concentration  is  related  to  the  logarithm  of  the  output  voltage,  it  was 
neccssar>'  to  calculate  the  concentration  for  each  sample  and  then  average 
the  concentrations  over  the  period  of  record  to  determine  the  true  mean 
concent  rat i on. 

Tlie  concentration  measurements  were  usually  started  as  near  the  bound- 
ary as  possible,  about  0.5  centimeter  above  the  crest  of  the  roughness 
elements.  nie  flume  was  stopped,  the  elevation  of  the  optical  equipment 
was  raised  and  recorded,  and  new  measurements  were  taken.  Tliis  procedure 
was  followed  until  an  elevation  was  reached  at  which  the  concentration 
was  too  low  for  the  optical  equipment  to  measure.  Tlie  optical  equipment 
was  then  lowered  in  a stepwise  fashion  until  near  the  boundary  to  obtain 
concentration  measurements  at  intermediate  elevations.  In  this  manner, 

7 to  15  concentration  measurements  were  obtained  to  describe  the  concen- 
tration distribution  for  one  flow  condition. 

3.  Resul ts . 

Sixty- five  concentration  distribution  curves  were  oiitained  of  the 
form , 

C = Cq  exp(M  Y)  , (9) 

where 

M = slope  of  the  curve  (in  feet“^) 

= concentration  at  the  base  elevation  (in  grams  per  liter") 

C = concentration  of  sediment  (in  grams  per  liter);  average  of 
tlie  concentrations  calculated  from  equation  (8)  for  eadi 
sample  of  record  during  the  period 

Y = elevation  aliove  the  crest  of  the  artificial  roughness  (in  feet) 
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'Hio  vai'iablo  used  to  dfscrilif  the  flow  eoiulition  is  the  flow 
\eloeity,  U , defined  as: 

IJ^,  = (.)  I.)/T  , flU) 

where  1.  is  the  amplitude  (in  feet)  and  is  equal  to  one-half  the  arc 
lengtti  which  j'asses  a fixed  point  during  one-half  cycle,  and  'I'  is  the 
period  of  oscillation  (in  secondsl.  'Hie  flow  conditions  studied  ranged 
from  a minimum  of  11^  = 0.235  foot  per  second  (minimum  to  allow  suspension 
of  sediment)  to  a maximum  of  = 1.18  feet  per  second  (maximum  allowed 
by  flume  construction).  The  amplitudes  and  periods  used  in  these  experi- 
ments ranged  from  0.235  to  1.60  feet  and  1.65  to  15.16  seconds,  respec- 
tively. Figure  5 shows  some  tvq'ical  concentration  distribution  curves 
tliat  were  obtained. 

For  the  65  different  flow  conditions  studied,  the  base  concentration, 
(eq.  y) , could  not  be  correlated  to  any  of  the  hydraulic  parameters 
but  depended  on  the  amount  of  sediment  in  the  flume.  For  different  flow 
conditions  the  sediment  in  the  flume  would  be  distributed  differently 
along  the  bottom  of  the  flume,  thereby  giving  a different  and  uncontrol- 
lable base  concentration. 

llie  65  experiments  also  showed  that  two  different  laws  exist  in  two 
ranges  of  conditions  governing  sediment  suspension.  For  am(ilitudes  of 
0.()93  foot  and  larger,  tlie  concentration  distribution  is  determined  by 
the  flow  velocity  alone.  For  amplitudes  less  tlum  0.693  foot,  the  con- 
centration distribution  is  a function  of  the  amplitude  relative  to  the 
wavelength  of  the  artificial  roughness. 

lliirty-six  of  the  65  concentration  distribution  curves  were  deter- 
mined asing  amiilitudes  of  0.693,  0.770,  0.925,  1.25,  and  1.60  feet.  For 
these  five  amplitudes,  it  was  found  that  the  slope  of  the  distribution 
curve  is  a function  of  the  flow  velocity,  independent  of  amplitude. 

Figure  6 graphically  illastrates  the  relationship  between  M,  the  slope 
of  the  concentration  distribution  curve  of  equation  (9) , and  U^,  the 
variable  of  equation  (10)  used  to  describe  the  flow  conditions  of  the 
flume.  These  data  arc  also  tabulated  in  Table  1 which  gives  the  periods, 
amplitudes,  and  use  of  the  optical  equipment.  llie  least  squares,  best 
fit  equation  for  this  relationship  is: 

M = -18.45  + (11.53)  % , (11) 

wliero  is  in  feet  per  second,  and  M is  the  slope  of  the  concentration 

distril'ution  curve  (in  fect”^).  lliore  is  no  statistical  evidence  to  indi- 
cate that  this  relationship  is  significantly  different  from  a higher  order 
polynomial . 

Results  of  experiments  by  Sliinoliara,  et  al.  (1958)  confirm  tlie  al'ove 
results.  Ihcy  found  tlie  same  linear  relationship  between  the  logarithm  of 
concentration  and  elevation  and  qualitatively  detennined  that  as  tlie  inten- 
sity of  the  flow  increased  the  slope  of  the  concent  nition  distribution 
curve,  M,  became  flatter. 
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Table  1.  (Concentration  distribution  data  for  \'g  = O.OoS  foot  per  second, 
amplitude  > 0.6'd3  foot. 


(Curve  n^Period  lAmnlitude 


518.01 

518.02 
518.0.5 

518.04 

518.05 


5 . 00 
7.34 
10.63 
12.70 
15.16 


IJ(5  Slope,  Variance  in 
d(ln(C)/d  least  squares 
(ft/s)  (ff^)  curve 


1.00 
0.681 
0.4  70 
0.394 
0 . 330 


-5.60 
-7.99 
-11. 28 
-13.91 
-15.39 


0.0182 
0.0184 
0.01 96 
0.02  3() 
0.0136 


Roma rks 


o;n  i cs 
stat i onarv 


602.01 

602.02 

602.03 

602.04 
602  .05 
602.06 
602.07 


612 

01 

612 

02 

612 

03 

612 

04 

612 

05 

705 

01 

705 

02 

705 

03 

705 

04 

707 

01 

707 

02 

707 

03 

707 

04 

8.8 
10.61 
12  .63 
15.10 

15.  14 
12.70 
10.68 
8.88 
6.15 


715.01 

713.02 

713.03 


925 
925 
0.925 
0.925 
0.925 
0.925 
0.925 


0.770 

0.770 

0.770 


1 . 

, 157 

0, 

, 737 

0. 

,504 

0, 

.420 

0, 

.348 

0, 

,292 

0, 

,2.35 

0, 

,•122 

0, 

,504 

0. 

,600 

0 

.721 

1 , 

.041 

1 , 

. 16 

0, 

,802 

0, 

,560 

0, 

.393 

0, 

. 968 

0, 

.671 

0. 

,4  70 

0 

.324 

-5 

96 

-12 

01 

- 1 1 

4 3 

-12 

90 

-14 

15 

-15 

60 

-14 

67 

1.18 

0.804 

0.558 


-7.05 

-8.88 

-13.12 

-14.20 


-5.89 
-8.58 
- 1 3. 65 


0.0628 

0.0140 

0.0151 

0.0043 


0.010 

0.034  1 


optics 
stat i onarv 


0, 

,00386 

0 

. 0333 

0. 

.0132 

0, 

.00359 

0.01 89 
0.00094 
0.00120 
0.00(1 89 


0, 

,00798 

0, 

.00194 

0 

.0(1,333 

0 

.00365 

0.0422 
0.00162 
0 . 00  1 1 3 


opt i cs 
stat i onarv 


opt i cs 
movi  ntj 


opt  i cs 
moving 


opt  1 cs 
moving 


717.01 

717.02 

717.03 

717.04 


0.770 

0.770 

0.770 

0.770 


0.392 
0 . 32  8 
0.468 
0.667 


•13.95 

-15.38 

■12.43 

-13.50 


0.0 -Sb 
0.00946 
0.0360 
0.0202 


op  tics 
mov ing 


72 1 .01 

721.02 

721.03 

721.04 


0.693 

0.693 

0.693 

0.693 


1.06 
0.72  3 
0.502 
0.352 


-6.29 

-9.40 

-13.10 

-16.30 


0.00913 

0.0023(1 

0.00514 

0.00320 


opt  i cs 
movi ng 


'Not  calculated. 


The  same  results  wore  obtained  more  recently  by  Kennedy  and  I.Qcher 
(1'.I72).  Their  investigation  examined  the  behavior  of  the  mean  sediment 
concentration  and  tl\c  periodic  sediment  concentration  fluctuations.  Ihe 
o.xperiments  wore  conducted  in  a st.ationary  flume  witli  a fixed  bed  on  which 
a limited  amount  of  loose  sedinx'nt  was  distributed.  Turbulence  for  sedi- 
ment siLspension  was  caused  by  surface  waves  generated  by  a wave  generator. 
Tlie  sediment  concentrations  wore  measured  with  optical  equipment  which 
incorporates  the  same  theoretical  principles  as  the  equij)ment  used  in  tiie 
swing-flume  e.xperiments  but  of  a mudi  smaller  site.  T!ic  smaller  size  and 
configuration  of  the  etiuipment  allowed  measurements  very  near  the  bed  and 
sampled  a much  smaller  volume  of  flow. 

Kennedy  and  I.ocher's  (1972)  experiments  in  mean  sediment  concentra- 
tion distribution  were  limited  to  a wave  height  of  0.24  foot,  a wave 
period  of  1.0  second,  and  a mean  water  depth  of  0.82  foot.  llie  mean 
sediment  concentration  was  measured  at  various  elevations  along  five 
evenly  spaced  verticals  in  the  flow.  The  spacing  of  the  verticals  was 
selected  to  cover  one  wavelength  of  the  bed  dune  shape.  A total  of  78 
data  points  was  measured,  and  when  the  logarithm  of  the  mean  concentra- 
tion was  plotted  against  elevation  above  the  bed  a well-defined  linear 
relationship  of  C = ex]’)(-36.5  Y)  was  obtained.  This  relationship  is 

identical  to  tliat  obtained  in  the  swing- flume  e.xperiments,  with  the  ex- 
ception of  the  high  rate  of  decay  of  sediment  concentration,  (-36.5). 
Kennedy  and  [.ocher  used  a quartz  sediment  of  0.  14-mi  1 1 imeter  mean  dia- 
meter in  their  experiments.  Tlic  settling  velocity  of  this  sediment  was 
not  reported;  it  was  probidily  about  0.050  foot  per  second,  which  is  43 
percent  greater  than  that  of  the  principal  plastic  sediment  (settling 
velocity  of  0.035  foot  per  second)  used  in  the  swing- flume  experiments. 
Swing-flume  measurements,  using  sediments  of  different  settling  velocities 
are  discassed  later  in  this  section.  For  a given  flow  velocity  a higher 
rate  of  decay  of  concentration  is  ex]iocted,  as  settling  velocity  increases 

Much  of  the  analysis  of  data  by  Kennedy  and  l.ocher  dealt  with  tlie 
periodic  sediment  concentration  fluctuations.  .Although  no  specific  con- 
clusions were  obtained  regarding  these  fluctuations,  the  data  did  indi- 
cate the  fluctuations  were  only  apparent  near  the  bed  (within  about  0.05 
foot  of  the  bed).  Ihis  exjilains  the  lack  of  periodicity  in  concentration 
fluctuations  for  the  swing-flume  measurements,  in  that  0.05  foot  is  near 
the  lower  limit  where  the  much  larger  optical  equipment  of  tlie  swing 
flume  could  be  used. 

Similar  exponential  concentration  curves  were  obtained  for  the  two 
different  methods  used  in  simulating  sediment  suspension  in  an  oscillat- 
ing flow;  i.e.,  an  oscillating  flow  over  a fixed  bed  as  used  by  Shinohara, 
et  al.  ( 1958)  and  Kennedy  and  [.ocher  (1972),  and  an  oscillating  l>ed  under 
a "stationary"  body  of  water  as  used  in  the  swing-flume  experiments. 

Figure  7 shows  the  relat  ionsiiip  between  M and  Uj?  for  two  condi- 
tions. One  set  of  tlata  represents  tlie  condition  that  the  o]iticai  equi['- 
ment  is  stationary  in  space;  the  other  <lata  [Joints  are  for  tlie  optical 
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equipment  moviuj;  with  the  flume.  All  otlicr  variables  in  these  cxiierimcnt' 
were  tlie  same.  ITic  amplitude  of  oscillation  was  U.'JdS  foot.  iliere  is  no 
statistically  significant  difference  in  the  results,  indicating  uniform 
su.spension  along  the  flume.  None  of  the  measurements  for  large  amjilitudes 
^greater  than  0.925  foot)  CTable  1)  wore  made  witli  tiie  optical  equipment 
moving  witli  the  flume,  llie  reason  for  this  is  at  large  amplitudes  the 
optical  equipment  would  move  into  an  area  of  the  fluid  where  secondary 
circulations  due  to  the  asymmetric  roughness  elements  may  exist  and  tiiere- 
fore  not  give  a representative  concentration. 

Tlie  remaining  29  of  the  65  ex]ieriments  were  made  using  amplitudes  of 
0.255,  0.465,  and  0.617  foot.  Table  2 is  a tabulation  of  the  data,  and 
Figure  S shows  the  results  of  the  65  experiments.  As  shown  in  Figure  8, 
there  is  a great  deal  of  scatter  in  the  data  for  small  amplitudes;  there- 
fore, only  qualitative  conclusions  have  been  made.  In  general,  the 
smaller  the  amplitude  the  smaller  the  slope  of  the  concentration  distri- 
bution curve 

l.xper iments , movies,  and  photos  were  used  to  determine  why  the  small 
amplitudes  do  not  obey  tlie  flow  velocit)'  relationship  of  the  larger  amf)li- 
tudes  given  in  ecpiation  (11).  Based  primarily  on  visual  observation,  the 
following  explanation  is  liypotlies  i zed . For  simill  amplitudes,  the  distance 
of  travel  of  tlie  flume  bottom  during  a half  cycle  is  not  great  enough  for 
the  boundary  layer  to  fully  develop  during  each  stroke  and  to  become  tur- 
bulent. Only  at  the  end  of  the  half  cycle  when  the  acceleration  forces 
cause  separation  is  sediinent  thrown  into  suspension.  Separation  only 
occurs  at  the  downstream  face  of  tlie  artificial  dunes.  Tlie  observed  sus- 
pension pattern  when  a 0.235- foot  amiilitude  of  oscillation  was  used  con- 
sisted of  plumes  of  suspension  separated  by  areas  of  zero  concentration, 
lltese  plumes  were  accentuated  because,  on  the  return,  tlie  half-cycle 
separation  at  the  downstream  face  occurred  such  that  the  succeeding  burst 
of  sediment  was  tlirown  into  ;ipproximately  the  same  region  of  fluid  as  in 
the  first  half  c>’cle;  i.e.,  the  amplitude  of  motion  was  about  equal  to  a 
multiple  of  the  wavelength  of  the  dune  shape.  iVlien  the  araiilitude  was 
increased  to  0.4(>5  foot,  a more  fully  developed  boundary  layer  was 
attained.  In  this  case,  separation  occurred  over  a somewhat  longer 
distance  of  travel  but  still  less  tlian  tiic  wavelength  of  the  dune  shape 
because  less  deceleration  force  was  required.  Hie  suspension  jiattern 
was  the  same  as  with  the  0.235-foot  amplitude  liut  much  less  distinct; 
the  plume.s  were  wider  and  overlapping.  Finally,  with  the  0.693-foot 
amplitude,  separation  occurred  over  a distance  equal  to  or  greater  tlian 
tlie  wavelength  of  the  dunes,  and  a uniform  longitudinal  concentration 
was  attained. 

Sfiecial  experiments  were  conducted  to  verify  the  above  hyjiotlies  i s . 

The  vertical  distribution  of  concentration  at  various  horizontal  loc.i- 
tions  in  the  fluid  was  measured  to  determine  if  the  concentration  distri- 
bution varied.  If  the  above  hypothesis  is  correct,  the  concentration 
distribution  shouhl  vary  along  the  horizontal  and  in  a regular  itumner 
determined  by  the  sliape  of  the  bed  dunes.  In  addition,  if  suspension 
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able  2.  (Concent  rat  i or  distribution  data  for  = 0.033  foot  per  second, 
amplitude  < 0.603  foot. 


Curve  Period  \nplitude 


621.01 
62 1 . 02 
621.05 

621.04 

62 1 . 05 


0.617 
0 . o 1 7 
0.617 
0.617 
0.1.17 


0 Slope  Variance  in 

dllnC)/d  least  squares  Remarks 
(ft/s)  (ff^)  curve 


1.170 
0.:’77 
0.557 
0 . 375 
0.263 


-11.15 
- 1(..  15 
-13.  10 
-15.12 
-10.05 


0 . 05 1 1 
0.0450 
0.01  34 
0.00268 
0.0233 


optl cs 
stat i onarv 


628, 

,01 

628, 

.02 

628, 

,03 

620 

.01 

620 

.02 

620, 

.03 

620 

.04 

711, 

,01 

71  1 , 

,02 

711, 

,03 

71  1 , 

,04 

711, 

,05 

711  , 

,06 

71  1 , 

,07 

802 

.01 

802 

.02 

802 

.03 

802 

,04 

802 

.05 

0.617 

0.617 

0.617 


814.01 

814.02 

814.03 

814.04 

814.05 


0.235 

0.235 

0.235 

0.235 

0.235 


0.465 
0.465 
0 . 465 
0.465 
0.465 


0.263 
0 . 376 
0.538 


0 . i’80 
1.177 
0.047 
0.448 


0.405 
0.585 
0 . 800 


0.570 
0.452 
0 . 363 
0.247 
0.207 


0.585 

0.585 

0.585 

0.585 

0.585 


-13.20 

-11.50 

-0.50 

-15.53 


0.16 
3.21 
6.63 
0 . 33 
-23.27 
-10.00 
-30 . 10 


-45.35 

-20.60 

-31.05 

-43.50 

-45.10 


0.00522 

0.0024^ 

0.0.305 

opt i cs 
moving 

0.00778 

0.00776 

0.00753 

0.00858 

optics 

moving 

0.0247 

0.00401 

0.00631 

0.00885 

0.0157 

0.00766 

0.0262 

opt i cs 
moving 

0.0150 

0.0204 

0.0292 

0.00577 

0.00821 

opt  i cs 
moving 

0 

,00481 

0 

.0239 

0 

.0208 

0 

.0152 

0, 

,0447 

optl cs 
station  a r\’ 
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for  small  amplitudes  is  a function  of  location,  the  shape  of  the  concen- 
tration distribution  curve  for  a particular  vertical  may  not  lie  the  same 
flume  velocity,  , hut  with  a large  amplitude.  Ilie  following  is  a 
description  of  the  experiments. 

Tlie  flume  was  adjusted  to  a 0.4(i5-foot  amplitude  and  a 3.  18-second 
period  was  selected.  A special  brace  with  attached  scale  was  construc- 
ted to  hold  the  optical  equipment  stationary  in  space,  and  therefore 
stationary  relative  to  the  fluid.  'Iliis  brace  allowed  relocation  of  tlie 
optical  equipment  to  any  desired  new  location  along  its  arc  while  main- 
taining tlie  freedom  to  change  the  elevation  of  the  optical  equipment 
relative  to  the  bed  of  the  flume.  An  initial  location  of  the  optical 
equipment  was  chosen  and  a vertical  concentration  distribution  measured, 
nie  optical  equipment  was  then  relocated  to  a new  position  and  another 
concentration  distribution  was  measured.  'Hiis  procedure  was  repeated 
until  five  concentration  distribution  curves  were  obtained,  l.adi  was  at 
a different  horizontal  location  relative  to  the  position,  at  a fixed 
phase,  of  the  flume  bed,  but  with  identical  flow  conditions  of  the  flume. 
Hie  horizontal  range  of  the  five  distribution  curves  was  about  equal  to 
one  wavelength  of  the  bottom  roughness  shape,  'llie  results  of  those  ex- 
periments (Fig.  9)  indicate  that  the  distribution  of  concentration  is  not 
uniform  along  the  flume  bed  for  the  0. 465-foot  amplitude  studied.  Figure 
9 also  indicates  that  the  rate  of  sediment  concentration  decay  is  less  in 
the  areas  between  the  dunes  and  greatest  at  the  crest  of  the  dunes.  Ibis 
irajilies  that  the  plumes  of  suspension  are  above  the  trough  of  the  bed 
shape,  where  it  would  be  expected  if  the  sediment  is  thrown  upward  and 
downstream  from  the  downstream  face  of  the  dimes. 

Figures  10  to  14  are  the  concentration  distribution  curves  measured 
for  the  locations  indicated  in  Figure  9.  Although  not  conclusive,  these 
data  indicate  that  the  vertical  distribution  of  sediment  in  the  plumes 
does  not  conform  to  the  exponential  distribution  of  equation  (9).  ibe 
relationships  shown  in  Figures  11  to  14  display  a slight,  hut  statisti- 
cally significant  curvature  and  were  obtained  at  locations  in  the  plume; 
whereas.  Figure  10  sliows  no  curvature  and  was  obtained  at  a location 
between  plumes.  Ibis  could  be  explained  by  the  existence  of  periodic 
stationary  eddies,  one  located  on  each  side  of  the  plume  axis.  Tlie  eddies 
would  sweep  sediment  into  the  base  of  the  plume  axis  giving  a relatively 
high  concentration  and  carry  sediment  out  of  the  top  of  the  plume,  causing 
the  concentration  to  bo  lower  than  would  be  predicted  in  a randomly  tur- 
bulent flow  field. 

Ibe  small  amplitudes  studied  represented  rare  prototype  conditions; 
therefore,  all  subsequent  experiments  were  limited  to  amplitudes  greater 
than  0.617  foot. 

4 . Fxperiments  Using  .Sediments  of  Different  Settling  Ve  loc  i t i es  . 

Two  different  sediments  of  the  same  black  plastic  were  selected 
and  investigated  to  determine  the  effect  of  settling  velocity  on  the 
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flume  is  reversing  its  direction  of  motion 


Figure  9.  M versus  location  of  optical  equipment  for  identical 
flow  conditions,  = 0.035  foot  per  second. 
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Y,  ELEVATION  ABOVE  CREST  OF  ROUGHNESS  ELEMENTS  (ft) 

I-'igure  10.  Concentration  distribution  for  run  814.01. 
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concentration  distribution.  'Flie  first  of  these  sediments  was  the  material 
wliich  passed  the  0.701-millimeter  sieve  and  was  retained  on  the  0.589- 
millineter  sieve.  Settling-velocity  measurements  of  225  randomly  cliosen 
particles  determined  the  moan  settling  velocity,  the  velocity  range,  and 
the  standard  deviation  as  0.0620,  0.0387  to  0.113,  and  0.0131  foot  per 
second,  respectively,  llie  second  sediment  was  the  material  wliich  passed 
the  0. 5S9-mi 1 1 imeter  sieve  and  was  retained  on  the  0.495-millimeter  sieve. 
Sett  1 ing- velocity  measurements  of  240  particles  of  this  material  deter- 
mined tlie  mean  settling  velocity,  the  velocity  range,  and  the  standard 
deviation  to  be  0.0498,  0.0285  to  0.0754,  and  0.00885  foot  per  second, 

< respectively. 

Table  3 gives  the  experimental  results  obtained  for  the  two  sediment 
types  and  Figure  15  shows  how  these  results  compare  with  the  results  for 
sediment  with  a settling  velocity  of  0.035  foot  per  second.  As  expected, 

L for  sediment  with  a higher  Vg,  the  concentration  of  sediment  decreases 

, with  elevation  above  the  bed  at  a higher  rate.  Only  this  qualitative 

conclusion  was  obtained.  Mot  enough  data  were  obtained  to  define  cpianti- 
F tativcly  the  relationship  between  Vg  and  '!. 

t Table  3.  Concentration  distribution  data  for 

Vg  = 0.0626  and  0.0498  foot  per  second. 


Curve 

No. 

Period 

(s) 

Ampi i tude 
(ft) 

Ur 

(ft/s) 

S 1 ope , 
d(  InC) /d 

(ft-i) 

Variance  in 
least  squares 
curve 

Vg  = 0 

.0626  ft/s 

913.01 

4.62 

1.25 

1.082 

-6.02 

0.0327 

913.02 

6.1)0 

1.25 

0.  758 

-11.90 

0.0154 

913.03 

9.43 

1.25 

0.530 

-15.68 

0.00497 

913.04 

11.46 

_ 1 

0.437 

-16.72 

0.000966 

o 

.0498  ft/s 

927.01 

3.  14 

0.848 

1.08 

-6.96 

0.00095 

927.02 

3 . 99 

0.848 

0.850 

-9.77 

0.00698 

927.03 

4.98 

0.848 

0.()82 

- 11.90 

0.00378 

927.04 

6.08 

0.848 

0.558 

-12.80 

0.00177 

llie  interesting  indication  of  these  measurements  is  that  the  increase 
in  the  rate  of  decay  with  elevation  of  the  concentration  for  a sediment 
of  higher  settling  velocity  is  not  as  groat  as  predicted  from  the  O'Brien 
(1933)  equation  for  concentration  equilibrium  conditions.  Consider  tlie 
concentration  equilibrium  equation  used  in  unidirectional  flow  for  sus- 
pended sediment: 

C Vg  + i;(dC/dY)  = 0 , (12) 


r-'igure  15.  M versus  \L  for 


1 

i 


where 


C = concentration 

V,  = settling  velocity  of  the  sediment  particle 

1;  = sediment  exchange  coefficient  Iwiiich  here  is  assumed  equal  to 

the  momentum  exchange  coefficient  in  unidirectional  flow) 

Y = elevation  above  the  bed 

Inserting  the  expressions  for  C and  (dC/dY)  obtained  from  equation 
(9)  into  equation  (12)  and  solving  for  I;  yield; 

f = -V'  /M  . (13) 

In  unidirectional  flow  and  for  relatively  low  sediment  concentrations, 
it  is  assumed  that  the  sediment  excliange  coefficient  is  equal  to  tlie 
momentum  exchange  coefficient,  is  a function  of  the  fluid  motion  only, 
and  is  independent  of  the  particle-settling  velocity.  For  this  assuniji- 
tion  to  be  valid  in  oscillating  flow,  the  slope  of  the  concentration  dis- 
tribution curve  must  be  directly  proportional  to  the  settling  velocity. 
Although  Figure  15  indicates  M is  proportional  to  V). , it  does  not 
indicate  direct  proportionality.  This  discrepancy  is  discussed  in 
Section  IV. 

5.  Summary  of  nx|3crimental  Results. 

Tlie  following  is  a summary  of  the  results  of  the  concentration  meas- 
urements and  a brief  iliscussion  of  their  limitations. 

a.  n>e  vertical  distribution  of  sediment  concentration  can  be  ex- 
pressed by  equation  (9).  Sediment  is  held  in  suspension  by  the  r.'uulom 
motion  of  turbulence  which  is  generated  in  the  boundary  layer  and  is 
transported  by  diffusioti  upward  while  decaying  continuously  because  of 
viscosity.  As  shown  later  in  Section  III,  the  turbulence  intensity 
decays  rapidly  with  elevation  above  the  bed.  Tlie  upper  elevation  to 
whicli  the  turbulence  can  diffuse  (the  water  surface  in  tlic  prototype 
and  the  wave  suppressent  board  in  the  flume)  is  larger  than  the  eleva- 
tion at  which  the  turbulence  intensity  decays  to  extremely  small  values; 
therefore,  this  upper  boundary  can  be  approximated  as  being  at  infinity. 
Under  these  conditions  it  is  rcasonalile  that  the  empirical  results  given 
by  equation  (9)  indicate  an  exponential  decay  (as  is  commonly  found  in  a 
diffusion  process),  and  only  become  zero  at  Y = infinity.  As  expected, 
the  concentration  distribution  in  oscillating  flow  is  different  than  in 
unidirectional  flow.  In  unidirectional  flow  tlie  turbulence  is  distrib- 
uted between  the  bed  and  the  water  surface,  the  turbulence  intensity  is 
significant  at  tlie  water  surface.  Because  the  water  surface  is  not  effec- 
tively at  infinity,  the  distribution  of  both  the  turbulence  and  the  sedi- 
ment concentration  would  be  different  than  in  oscillating  flow. 
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b.  Hie  base  concentration,  feq.  y),  for  flume  measurements 

is  a function  of  the  sediment  charge  in  the  flume  and  therefore  could 
not  be  correlated  to  flow  hydraulics.  As  shown  by  liinstein  [1950),  a 
flow  is  only  capable  of  transporting  a limited  amount  of  sediment  of  a 
given  sicc.  Tliis  limiting  capacit)’  is  determined  by  the  flow  velocity. 
Sediment  characteristics,  and  roughness  of  the  boundary.  In  addition, 
the  flow  will  only  transport  this  capacity  rate  if  there  is  a sufficient 
supply  of  sediment  available.  Otherwise,  tlie  transport  rate  will  he 
reduced  by  tiie  ratio  of  the  supply  rate  to  the  capacity  rate.  Because 
the  capacity  transport  rate  for  a given  flow  is  determined  by  the  proba- 
bility of  a particular  sediment  particle  being  subjected  to  sufficient 
hydraulic  forces  to  move  it,  there  must  be  some  particles  in  the  bed  that 
are  not  in  motion  at  any  instant  of  time.  Had  there  been  enough  sediment 
in  the  flume  to  satisfy  the  flow's  sediment  transport  capacity,  the 
measured  Cq  could  have  been  correlated  to  flow  velocity.  Unfortunately, 
under  those  conditions  some  sediment  must  be  loosely  deposited  on  the 
flume  bed,  theretiy  changing  the  fixed-bed  geometry  and  roughness.  Tliere- 
fore,  only  the  flow's  capacity  to  transport  sediment  of  a specific  size 
can  be  estimated.  This  estimate  must  use  the  capacity  base  concentration 
calculated  from  Kalkanis'  (1964)  bedload  equation  and  not  the  base  concen- 
trations measured  in  this  investigation. 

c.  In  tlie  range  of  flow  velocities,  0.2  foot  per  second  < < 1.1 

feet  per  second,  for  amplitudes  of  oscillation  equal  to  or  greater  than 
0.693  foot,  and  for  Vg  = 0.035  foot  per  second,  the  slope  of  the  exponen- 
tial distribution  of  sediment  concentration  is  a function  of  flow  veloc- 
ity only,  llic  slope,  M,  can  be  approximated  from  the  flow  velocity, 

Uq,  by  equation  (11).  Tliis  equation  is  only  a best  fit  empirical  rela- 
tionship and  cannot  give  reasonable  approximations  of  M for  Uq  values 
very  far  outside  the  stated  range.  Tliis  becomes  apparent  when  substitut- 
ing in  a large  value  of  e.g.,  2.0  feet  per  second,  and  calculating 

M.  The  result  would  be  a positive  value  for  M;  i.e.,  the  concentration 
of  sediment  increases  with  elevation  which  is  not  reasonable.  Tlie  limit- 
ing value  of  M for  extremely  large  values  of  should  bo  zero,  or 

uniform  concentration  of  sediment  throughout  the  depth.  At  the  other 
extreme,  equation  (11)  gives  a value  of  M = -18.45  feef^  for  = 0 foot 
per  second  which  is  also  not  reasonable.  But,  for  low  values  of  U^,  any 
continuous  function  is  not  expected  to  give  a correct  relationship  since 
at  some  point  the  flow  changes  from  turbulent  to  laminar.  In  laminar 
flow  there  is  no  turbulence  and  therefore  no  sediment  suspension.  As 

is  increased,  the  flow,  at  some  velocity,  suddenly  changes  from  laminar 
to  turbulent  and  just  as  suddenly  the  suspended-sediment  concentration 
changes  from  zero  to  some  positive  value.  Ttiereforc,  the  relationship 
exi^ressed  by  equation  (11)  becomes  invalid  at  some  small  value  of 

d.  For  the  flow  velocities  studied,  the  sedi ment -sett  1 i ng  velocity 
has  a significant  effect  on  the  slope  of  the  concentration  distribution 
curve.  Not  enougli  data  were  olitained  to  define  the  relationship  between 
Vg  and  M,  but  the  data  did  yield  the  qualitative  relationship  that  for 
constant  , M decreases  (or  becomes  a larger  negative  value)  with  in- 
creasing Vg . As  discussed  earlier,  if  the  sediment  exchange  coefficient 
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is  equal  to  the  momentum  exchange  coefficient,  as  is  assumed  for  unidirec- 
tional flow,  the  slope,  M,  should  be  directly  proportional  to  VV . The 
experimental  results  for  oscillating  flow  did  not  indicate  direct  propor- 
tionality. Therefore,  the  sediment  exchange  coefficient  for  oscillating 
fli'W  given  by  equation  (13)  is,  as  yet,  an  undetermined  function  of  Vg 
and  A discussion  in  Section  IV  indicates  why  this  result  is  expected. 

III.  DISTRIBUTIONS  OF  TURBULENT  VELOCITY  FLUCTUATIONS 
1.  l.xperimcntal  Apparatus. 


Successful  measurements  of  turbulent  velocity  fluctuations  in  fluids 
have  been  imide  using  either  constant  current  or  constant  temperature 
hot-film  anemometers.  A constant  temperature,  quartz-coated  hot-film 
sensor,  model  number  6010  made  by  Tltermo-Systems  Incorporated,  Minneapolis, 
Minnesota,  was  used  in  this  investigation.  The  sensor  was  connected  to  a 
1050  series  anemometer  also  made  by  Thermo-Systems  Incorj)orated.  Tlie 
anemometer  uses  a bridge  and  feedback  system  to  maintain  a constant  resis- 
tance and  therefore  a constant  temperature  of  the  sensor.  Any  change  that 
affects  the  heat  transfer  between  the  sensor  and  the  environment  is  reflec- 
ted in  tlie  voltage  output  of  the  bridge.  Tliis  output  voltage  is  amplified 
and  recorded  on  magnetic  tape.  Tlie  record  of  voltage  fluctuations  is  then 
converted  by  use  of  a calibration  curve  to  a record  of  velocity  fluctua- 
tions. A schematic  of  the  hot-film  bridge  is  shown  in  Figure  16.  Tlie 
hot-film  sensor  and  probe  are  shown  in  Figure  17. 

Calibration  of  the  hot-film  anemometer  was  done  in  the  calibration 
tank  (Fig.  18)  which  was  divided  into  two  chambers,  a fore  chamlicr  and  a 
calibration  chamber.  The  fore  chamlier  contained  two  wire  screens  to 
ensure  a uniform  velocity  distribution.  Tlie^  two  chambers  were  connected 
by  a 1-inch-diameter  nozzle  located  at  the  midpoint  of  tlie  partition 
separating  the  two  chambers,  5.25  inches  above  the  bottom  of  tlie  tank. 

An  overflow  was  located  at  the  downstream  end  of  the  calibration  cliamber, 
4.75  indies  above  the  nozzle.  Water  was  supplied  from  the  deaeration 
tank,  the  rate  controlled  by  a 0.25-inch  needle  valve.  'Hie  flow  into 
the  calibration  chamber  was  a submerged  jet.  Tlie  probe  with  sensor  was 
held  at  the  downstream  face  of  the  nozzle  by  clamps  connected  to  a point- 
gage  assembly.  Tlie  point-gage  assembly  was  used  to  raise  and  lower  the 
sensor  known  amounts  to  obtain  velocity  measurements  across  the  diiimeter 
of  the  nozzle.  'Itie  water  collected  at  the  overflow  in  a measured  time 
was  weighed  to  determine  the  flow  rate  and  moan  nozzle  jet  velocity. 

Measurements  of  the  turbulent  velocity  distributions  were  conducted 
in  the  swing  flume.  Only  minor  modifications  to  the  flume  apparatus 
were  necessary  to  accommodate  the  anemometer  equi]iment. 

As  shown  in  Figure  17,  the  sensor  is  extremely  delicate  and  there- 
fore cannot  be  used  in  flows  containing  solid  particles.  For  this  reason 
the  flume  had  to  be  cleaned  of  all  the  black  plastic  sedimiuit  used  in 
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BRIDGE  VOLTAGE 


Figure  Ifa.  Functional  'ichematic  of  iiot-wirc  liridge  circuit 


concentration  measurements.  To  ensure  that  the  natural  sediment  glued 
to  the  artificial  dunes  did  not  break  loose  and  damage  the  sensor,  the 
surface  of  the  dunes  was  sprayed  with  a tliin  film  of  plastic.  Tlie  plas- 
tic was  thin  enough  to  not  alter  the  surface  roughness  and  strong  enough 
to  hold  the  sediment  in  place. 

To  accommodate  the  anemometer  equipment,  the  yoke  support  used  in 
concentration  measurements  had  to  be  removed  from  the  support  frame 
(Fig.  3)  . Ttie  liot-film  probe  assembly  (Fig.  19)  was  attached  to  the 
support  frame  at  a position  midway  across  the  flume.  Ttie  elevation  of 
the  sensor  could  be  changed  by  loosening  two  friction  clamps  and  repo, i- 
tioning  the  probe  holder  in  relation  to  a scale  fixed  to  the  assembly. 

A 6-inch-square  opening  was  cut  into  the  wave  suppressent  board  to  allow 
the  sensor  to  be  lowered  to  the  flume  bottom. 


experimental  Procedure. 


riie  hot-film  sensor  is  extremely  sensitive  to  both  water  temj^erature 
and  water  quality.  Ttierefore,  it  was  necessary  to  determine  a new  cali- 
bration curve  each  day  that  velocity  measurements  were  made  in  the  flume. 
Water  temperature  readings  made  during  calibration  measurements  were  com- 
pared to  readings  made  during  flume  measurements  and  did  not  vary.  It 
was  necessary  to  use  deaerated  water  in  both  the  calibration  and  flume 
measurements  to  prevent  air  bubbles  from  adliering  to  the  hot- film  sensor, 
thereby  either  burning  out  the  sensor  or  altering  its  lieat  transfer  char- 
acteristics. Periodically,  during  both  calibration  and  flume  measurements 
a record  was  made  of  the  base  voltage;  i.e.,  the  voltage  for  the  sensor  in 
still  water.  Ttiese  voltage  readings  were  then  averaged  to  obtain  a mean 
base  voltage  which  is  needed  in  calculating  the  calibration  curve. 


Both  days  that  velocity  measurements  were  made  with  the  anemometer, 
a 5-percent  o.verheat  of  the  sensor  was  used.  Tliis  ensured  uniform  sensi- 
tivity of  the  sensor  and  base  voltage  of  the  measurements.  Tlie  conditions 
for  measurements  on  both  days  were  so  similar  that  the' two  calibration 
curves  could  not  be  distinguished.  This  allowed  all  the  calibration  data 
to  be  used  for  one  curve  and  only  one  equation  used  to  convert  voltage 
into  velocity. 


a.  Calibration  of  the  Hot-Film  Sensor.  Tlie  hot- film  sensor  was 
positioned  at  the  center  and  as  close  to  the  downstream  face  of  the  cali- 
bration nozzle  as  possible  (approximately  one-eighth  inch).  A liigh  flow- 
rate through  the  flume  was  obtained  by  fully  opening  the  needle  valve. 
This  flow  rate  was  allowed  to  continue  until  a steady  flow  through  tlie 
calibration  tank  was  est.iblished  (about  5 minutes),  'file  flow  through  the 
nozzle  was  monitored  by  the  continuous  voltage  output  of  the  hot-film 
bridge  displayed  on  the  anemometer  equipment.  IVlien  the  voltage  readings 
became  constant  with  time  the  flow  rate  was  steady.  A 2't-second  magnetic 
tape  record  was  then  made  at  the  rate  of  117  samples  per  second.  Simul- 
taneous with  the  voltage  record,  a flow  rate  measurement  was  made.  Tlie 
flow  rate  was  determined  by  timing  the  period  required  to  collect  10  to 
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15  pounds  ot-'  water  at  the  overflow  of  the  calibration  tank.  Tlie  water 
was  weighed  on  a scale  to  the  nearest  0.05  jjound.  'Die  flow  through  the 
flme  was  then  reduced  by  partially  closing  the  needle  valve  and  a second 
set  of  measurements  made.  Tlie  procedure  was  repeated  until  five  voltage- 
velocity  points  were  measured  for  a calibration  curve. 

'I'he  results  of  the  two  calibration  curves  are  shown  in  Figure  20. 

■file  velocity  range  measured  was  from  0.0374  to  0.580  foot  per  second. 

'Ilie  relationship  between  velocity  and  voltage  was: 

ln^(F'  - ho)  = (0.422)  ln,;(Ug)  + 7.147  , (14) 

where  li'  is  proportional  to  tlie  output  voltage  of  the  hot-film  bridge, 
and  is  pro|)ortional  to  the  mean  output  voltage  for  the  sensor  in 

still  water.  U , the  actual  velocity  of  the  fluid,  was  calculated  by 
dividing  the  mean  velocity  determined  from  flow  rate  measurements  by  the 
nozzle  coefficient,  C^,.  'Hie  log-log  relationship  shown  in  Figure  20  is 
consistent  with  the  theoretical  results  from  the  manufacturer  and  the 
results  by  Uas  (19o8).  Because  magnitudes  of  the  voltage  readings  for  a 
given  flow  velocity  are  dependent  on  the  water  temperature,  water  quality 
overheat  percentage,  and  amplification  of  the  sensor  output  signal,  no 
attempt  was  made  to  compare  quantitatively  the  measured  calibration  curve 
witli  other  published  curves. 


b.  Oeterminat  ion  of  the  Nozzle  Coefficient.  Tlie  sensor  was  posi- 
tioned close  to  the  downstream  face  and  at  the  lower  edge  of  tlie  nozzle. 
After  establishing  a high,  cons  ant  flow  rate  through  the  nozzle,  the 
voltage  from  tlie  liot-film  bridge  was  recorded.  Without  interrupting  the 
flow,  the  sensor  was  then  raised  a small  measured  amount  and  a second 
voltage  record  made,  lliis  procedure  was  repeated  until  the  sensor  was 
at  the  upper  edge  of  the  nozzle.  Tlie  velocity  of  flow  was  estimated  for 
each  location  from  a calibration  curve  similar  to  Figure  20  but  not  cor- 
rected by  a nozzle  coefficient.  Tlie  velocity  profile  across  the  diameter 
of  the  nozzle  for  the  high  flow  rate  is  shown  in  Figure  21(a).  Because 
of  the  high  flow  rate,  the  water  surface  elevation  in  the  deaeration  tank 
was  lowering  a significant  amount,  thereby  decreasing  the  flow  rate 
tlirough  the  caliliration  tank.  Tliis  is  the  reason  for  the  lower  measured 
velocity  at  tlie  toji  of  the  nozzle.  Tlie  mean  velocity  was  determined  by 
integrating  the  velocity  ]irofile  across  the  jet.  The  nozzle  coefficient 
was  then  calculated  from: 


(.'  = [(V  ) (dia.)^  ,]/[(dia.)2 

n mecPi  'jet'  'nc 


where 


Cr.  = nozzle  coefficient 


'^rne,T>-  ~ mean  velocity  across  the  jet 


(dia.)j^^  = diameter  of  the  jet 

(dia.)^j^,,  = diameter  of  the  nozzle  (1  inch) 

(V  ) = velocity  determined  from  voltage  measurements 

' neae . - ^ 
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Figure  21.  Velocity  profile  across  calibration  nozzle. 


riio  nozzle  coefficient  for  the  high  velocity  (U.569  foot  jier  secoiulj  uas 
tiien  calculated  to  he  0.99.  'Ihis  set  of  raeasurci.ients  was  repeated  for 
a low  flow  rate.  'Hie  results  of  these  measurements  arc  shown  in  figure 
dl(b).  Hie  nozzle  coefficient  for  the  low  vclocit)'  (0.072  foot  per  sec- 
ond) was  calculated  to  he  0.98.  llic  nozzle  coefficient  for  intermediate 
velocities  was  assumed  to  be  a linear  interpolation  between  the  two  altove 
values . 

c . Flume  Measurements  of  Turbulent  Velocity  Distributions.  'Hi e 
velocity,  which  is  of  significance  in  attempting  to  describe  the  sus- 
pension of  sediment  in  an  oscillating  flow,  is  the  vertical  component  of 
vclocit)'  fluctuations  due  to  turbulence,  v'.  Ulien  the  swing  flume  is 
0])crating,  the  only  velocities  which  exist  in  the  fluid  are;  (a)  tlie 
three  directional  components  of  velocity  fluctuations  caused  liy  turbu- 
lence, and  (b)  tlie  oscillating  flow  contained  in  the  boundary  layer. 

Tlie  boundary  layer  extends  only  a few  millimeters  above  the  bed  of  tlie 
flume.  Sediment  is  in  suspension  at  an  elevation  considerabl)'  above 
the  upper  limit  of  the  boundary  layer.  Tliercfore,  the  flow  regime  of 
interest  lias  no  measurable  mean  or  periodic  velocities,  only  the  random 
motions  Caused  by  the  turbulence  diffusing  upward  from  the  bed.  'Hie 
problem,  then,  is  to  determine  only  the  vertical  component  of  the  veloc- 
ity fluctuations. 

Das  (19()8)  developed  a method  of  measuring  v'  in  a still  Inxiy  of 
water  with  an  oscillating  rough  bed.  'Iliis  method  involved  imparting  an 
oscillating  motion  to  the  sensor  in  the  vortical  direction.  if  only 
measurements  made  during  the  peak  velocity  of  the  sensor  are  considered, 
then  the  total  velocity  affecting  heat  transfer  from  the  sensor  is: 

V|  = (V  + v')^  + u'2  + w'2  , ( !(,) 

where 

Vg  = total  velocity  affecting  heat  transfer 

V = peak  velocity  of  the  sensor;  known  from  the  period  and 

amplitude  of  osci 1 lat ion 

v'  = vertical  component  of  tlic  turbulent  velocity  fluctuations 

u'  and  w'  = turlnilent  velocity  fluctuations  in  the  remaining  two 
d i rect i ons 

Dividing  this  equation  by  \’^  )•  i e Ids  : 

(V2/V’2)  = I + (2  v'/V)  + (v'/V)"  + (u'/V)^  + (w'/\l-  • (n 

If  the  sensor  oscillation  is  such  that  \'  ■■  v',  u',  ami  w',  then  tiu' 
above  equation  can  be  approximated  by: 

fV,2/V?)  = 1 1-  (2  v'/V) 


( 1S| 
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lliis  equation  is  sensitive  to  only  v'.  Fixperiments  by  Das  were  conduc- 
ted in  a stationary  flume  in  wliich  only  a horizontal  bottom  plate  was 
oscillated  to  produce  turbulence.  Tlie  results  showed  some  promise  for 
the  method. 

Wlien  Das'  method  of  measuring  v'  was  tried  in  the  swing  fliwie,  it 
was  not  successful  because  of  excessive  vibrations  of  the  sensor.  Tliese 
vibrations  were  mainly  due  to:  (a)  the  long  holder  required  to  extend 

the  sensor  to  the  flume  bottom,  and  (b)  attaching  this  holder  to  the 
support  frame  whicli  was  indirectly  subjected  to  the  vibrations  from  the 
flume  motion.  Tlie  velocities  of  the  sensor  due  to  vibrations  were 
greater  than  the  velocities  of  the  turbulent  flow;  therefore,  no  com- 
ponent of  turbulent  velocity  could  be  distinguished. 

,\n  ai>proximat  ion  to  v'  had  to  be  obtained  based  on  the  following 
assumptions.  It  was  assumed  that  at  a given  elevation  the  root-mean- 
square  value  of  the  throe  components  of  turbulent  velocity  fluctuations 
is  proportional  to  each  other.  It  was  also  assumed  that  the  heat  con- 
vected  from  the  sensor  due  to  velocities  in  the  direction  parallel  to 
the  sensor  axis  was  insignificant  compared  to  the  heat  convected  by 
velocities  perpendicular  to  the  axis.  This  assumption  is  justified  in 
that  the  hot-film  has  directional  properties  making  the  maximum  sensi- 
tivity at  right  angles  to  the  flow.  Also,  the  aspect  ratio  (length- 
diameter)  of  the  sensor  is  such  that  its  properties  approach  those  of 
an  infinite  wire  where  there  is  no  effect  of  a longitudinal  velocity. 
Bascil  on  these  assumptions,  the  effective  velocity  causing  heat  convec- 
tion is,  as  an  average: 

Ilg  = [v'2  (K'2  + 1)]'^  , (19) 

where  is  the  velocity  corresponding  to  the  output  voltage  of  tlie 

hot-film  bridge,  k'  is  the  constant  of  proportionality  between  the 
vertical  component  and  one  of  the  horizontal  components  of  turbulent 
velocity  fluctuation,  and  v'  is  the  vertical  comjionent  of  the  turbulent 
velocity  causing  heat  convection  from  the  sensor.  It  is  apparent  then 
that  the  sensor  must  be  placed  in  the  flume  with  its  axis  horizontal. 

Ilie  magnitude  of  the  vertical  component  of  velocity  fluctuation  can  then 
be  calculated  and  is,  as  an  average: 

v'  = [u2/(K'2  + n r*  . (dO) 

Although  equation  (20)  is  only  an  approximation,  the  assumptions  used 
do  not  affect  the  basic  relationships  (a)  between  the  root -mean-S(|u:ire 
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value  of  v'  and  elevation  above  the  bed,  and  (b)  between  the  root- 
mean-square  value  of  v'  at  a fixed  elevation  and  the  flow  velocit\-,  U^. 
Tlie  assumptions  also  allow  an  approximation  of  the  absolute  magnitude  of 
the  root- mean-square  value  of  v'.  Tlie  horizontal  component  of  turbulent 
velocity  fluctuation  is  probably  on  the  same  order  of  imignitude  as  the 
vortical  component  and  therefore,  for  qualitative  analysis,  the  value 
of  K’  in  equation  (20)  can  be  approximated  as  equal  to  unity. 


Tile  procedures  used  in  measuring  velocity  distributions  in  the  flume 
wore  as  follows.  A period  and  an  amplitude  of  flume  oscillation  were 
selected  and  the  flume  linkage  adjusted  to  give  a symmetric  motion. 

After  the  flume  was  filled  with  deaerated  water  to  the  elevation  of  the 
wave  suppressent  board,  the  asymmetric  roughness  elements  were  adjusted 
to  eliminate  secondary  currents  in  the  central  part  of  the  flume.  Tlie 
sensor  was  then  placed  in  the  flume  as  near  the  bottom  as  possible  and 
its  elevation  recorded.  Tlie  flume  was  started  and  the  motion  allowed  to 
continue  until  equilibrium  flow  conditions  were  established.  A record 
of  the  hot-film  bridge  output  voltage  was  made  on  magnetic  tape,  tlie 
length  of  which  was  an  integer  multiple  of  the  flume  oscillation  period. 
Tlie  flume  was  stopped  and  the  sensor  elevation  raised  for  a new  measure- 
ment. Tlic  procedure  was  repeated  until  an  elevation  was  reached  at  which 
the  velocity  fluctuations  were  too  small  to  be  accurately  measured  with 
the  anemometer.  Tlie  sensor  was  then  lowered  in  a stepwise  manner  to 
obtain  velocity  measurements  at  intermediate  elevations.  In  tills  mnner, 
10  to  13  velocity-elevation  measurements  were  obtained  to  give  a velocity 
distribution  for  the  flow  condition  used.  The  period  of  the  flume  was 
changed  and  the  measurements  repeated  to  give  a second  velocity  distribu- 
tion. In  all,  four  velocity  distributions  were  obtained  for  four  differ- 
ent flow  conditions. 


3.  Results. 


The  purpose  of  the  velocity  measurements  was  to  obtain  the  following 
three  relationships  needed  for  an  analy'sis  of  the  suspended- load  e<iuation; 
(a)  An  approximation  of  the  magnitude  of  the  root-mean -square  value  of  v' 
versus  flow  velocity,  Uq;  (b)  the  distribution  of  tlie  root-mean-square 
value  of  v'  versus  elevation  above  the  bed;  and  (c)  the  distribution 
with  time  of  v'  at  a constant  location  in  space.  Results  pertaining 
to  the  third  unknown  listed  above  will  be  discussed  first. 

i\vo  sets  of  data  were  analyzed  to  determine  tiie  distribution  with 
time  of  v'.  Hie  period  and  amplitude  of  flume  oscillation  for  both  sets 
of  data  were  10.48  seconds  and  0.925  foot,  respectively.  In  botli  cases, 
the  length  of  record  analyzed  was  10.48  seconds  (1 , 224- vo 1 tage  samples). 
One  sot  of  data  was  taken  at  an  elevation  of  0.  I(i8  foot  aliove  the  crest 
of  the  artificial  dunes,  the  other  0.209  foot  above. 

Hie  data  were  analyzed  in  the  following  manner.  For  each  voltage 
samiilc  recorded,  the  effective  heat  transfer  velocity,  U , , was  calcu- 
lated from  equation  (14).  Ilie  velocities  were  ordered  and  percentages 
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equal  to  or  less  than  various  selected  velocities  calculated.  Tliese 
percentages  were  divided  by  two  and  plotted  against  velocity  on  normal 
probability  paper.  llie  percentages  were  divided  by  two  in  order  to 
atljust  for  the  fact  that  the  anemometer  measured  the  absolute  effective 
velocity  without  regard  to  its  direction;  i.e.,  the  velocities  in  each 
range  were  composed  of  an  equal  number  of  negative  and  positive  veloci- 
ties, thereby  giving  twice  the  percentage  of  actual  positive  velocities. 

As  shown  in  figure  22(a,  b) , those  plots  approximate  straight  lines  and 
tlie  50-percent  velocity  is  zero,  thereby  indicating  that  the  distribu- 
tion of  the  turbulent  velocity  fluctuations  is  approximately  normal  with 
a moan  of  zero.  Similar  data  by  Uas  (1968)  give  the  same  results. 

'Hie  above  result  suggests  that  the  standard  deviation,  s,  of  the 
normal  distribution  (which  equals  tiie  root-mean-squarc  velocity)  be  used 
as  tlie  velocity  scale  describing  turbulence  intensity  for  a given  eleva- 
tion. 

Tlie  data  were  then  analyzed  to  determine  the  distribution  of  s witli 
respect  to  elevation.  For  each  sample  of  a voltage  record,  the  effective 
heat  transfer  velocity  was  calculated  from  equation  (14).  From  the  velo- 
city record,  the  root-mean-square  effective  velocity  was  calculated. 
Knowing  this  velocity  the  velocity  scale  for  the  elevation  at  which  the 
record  was  made  was  calculated  from: 

s = UJ/vT  , (21) 

where  s is  the  velocity  scale  and  is  etiual  to  the  standard  deviation  of 
v'  for  tlic  elevation  of  the  record  and  flow  conditions  of  the  flume,  11^ 
is  the  root-mean-square  effective  heat  transfer  velocity  of  the  record, 
and  /J  comes  from  equation  (20)  when  K'  is  assumed  equal  to  unity. 

llie  velocity  scale  was  plotted  against  elevation  on  semi logar i thmi c 
paper  to  give  the  relationships  shown  in  Figures  23  to  26.  In  general, 
tliese  relationships  can  be  exjiressed  by: 

s = s^  exp (A  V)  , (22) 

where  S|5  is  tlie  value  of  the  velocity  scale  (in  feet  per  second)  at  the 
elevation  of  the  crest  of  the  artificial  bed  dunes,  A is  the  slope  of 
the  exponential  curve  (in  feet"^),  and  Y is  the  elevation  (in  feet) 
above  the  crest  of  the  bed  dunes.  Tlie  flume  flow  conditions  for  which 
a velocity  scale-elevation  distribution  was  measured  were  1)^  = 0. 3.5.3, 
0.510,  0.748,  and  0.930  foot  per  second. 

Comparison  of  the  four  velocity-elevation  distributions  revealed  that, 
for  the  range  of  flow  conditions  studied  and  bed  roughness  used,  the  slope 
of  the  cxiioncnt  i al  relationship  was  constant.  Tliis  implies  that  tlie  in- 
tensity of  the  turbulence  decreases  in  a manner  which  is  independent  of 
the  flow  velocity  generating  the  turbulence.  Hie  constant  rate  of  velo- 
city dec.'iy  appears  to  apjily  to  elevations  near  the  bed,  within  1.5  centi- 
meters as  neasureil  with  the  hot-film  sensor. 
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TO  OR  GREATER  THAN  VELOCITY  INDICATED 
(NORMAL  PROBABILITY  SCALE) 


The  base  velocity,  s^  (cti-  22),  was  found  to  be  a function  of  tiic 
flume  flow  conditions.  'Hie  relationsliip  between  s^  and  U , is  shown 
in  figure  27.  from  ti\e  limited  data  it  was  onl\-  possible  to  determine 
an  approximate  mathematical  relationship  lietween  s and  U^.  This 
relationship  is: 

s^  = (0.0885)  + 0.0557  , (23) 

which,  from  boundary  considerations,  only  gives  approximate  s^  values 
in  the  range  of  experimental  values  of  . 


Summarv  of  Hxperimental  Results. 


ilie  following  is  a summary  of  the  results  of  the  turbulent  velocity 
fluctuation  measurements  and  a brief  discussion  of  their  limitations. 


a.  ’llie  velocity  fluctuations  caused  by  turbulence  are,  for  a con- 
stant elevation  above  the  bed  and  a constant  flow  velocity,  approximately 
normally  distributed  with  a mean  of  zero.  The  standard  deviation  of  the 
distribution,  s,  is  used  as  the  velocity  scale  to  measure  turbulence 
intensit)’  at  any  elevation.  Most  turbulent  velocity  fluctuation  measure- 
ments give  approximately  Gaussian  results,  although  it  is  known  that 
except  for  isotropic  turbulence,  the  distribution  cannot  be  Gaussian. 

b.  for  the  range  of  flow  conditions  studied  and  the  bed  roughness 
used,  the  velocity  scale  can  be  expressed  by  equation  (22).  The  exponen- 
tial nature  of  tb.  is  relationship  was  discussed  in  Section  11,  paragraph 
4(a).  llie  relationship  also  conforms  to  tlie  boundary  conditions  of  the 
flow.  .\s  expected,  the  turbulence  intensity  assumes  a limiting  value, 

S|5 , at  the  ocean  l)ottom  (Y  = 0).  Tills  limiting  value  of  turbulence  in- 
tensity is  determined,  in  some  manner,  by  tlie  flow  velocity,  U^.  As 
the  turbulence  diffuses  upward  its  intensity  decays  because  of  viscosity. 
Tlie  body  of  fluid  into  which  the  turbulence  diffuses  is,  by  comparison, 
e.xtremel)’  large;  therefore,  tlie  empirical  relationship  is  expected  to 
indicate  that  the  turbulence  intensity  decays  to  zero  at  an  infinite 
distance  from  tlie  bed. 


c.  llie  slope.  A,  of  the  exponential  distribution  of  the  velocity 
scale  is  constant  with  respect  to  elevation  and  constant  throughout  the 
range  of  flow  conditions  studied.  It  was  found  to  be  -10.57  feet 

lliis  result  is  not  suiqirising  since  the  rate  of  turbulence  intensity 
decay  for  the  oscillating  flow  conditions  measured  is  determined  by  vis- 
cosity. Therefore,  for  fluids  of  the  same  viscosity  and  density,  the 
rate  of  decay  should  be  constant  and  independent  of  flow  velocity. 

d.  llie  base  velocity  scale,  s^  (e<|.  22),  is  a function  of  the 
flow  velocity  ami  can  be  approximated  by  equation  (23).  lliis  equation 
is  a best  fit  relationsliip  of  the  empirical  data  and  does  not  appl>’  for 
flow  velocities  outi.iile  the  measured  range.  This  becomes  obvious  by 
letting  11,5  = 0 foot  ]ier  second  and  finding  s = ti.0557  foot  per  second. 
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(ft/s) 


1.0 


Sq  - VERTICAL  COMPONENT  OF  TURBULENT  VELOCITY 
FLUCTUATIONS  AT  THE  CREST  OF  THE  BED  ROUGHNESS 
ELEMENTS  (ft/s) 

ligiin.'  J7.  Base  vertical  turlnilcnt  velocity  scale  versus  flume 
velocity  for  a constant  amiilitmle  = ().92ri  foot. 
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In  tlie  range  of  flow  velocities  where  equation  (23)  is  valid,  the  calcu- 
lated Sq  is  an  approximation  of  a base  volocit\'  at  an  arbitrary  eleva- 
tion, the  crest  of  the  bed  dunes.  To  apply  the  experimental  results  to 
a real  situation,  the  const;uits  in  equation  (23)  must  be  adjusted  to  give 
So  values  at  the  elevation  of  the  top  of  tlie  bedload  layer,  lliis  eleva- 
tion depends  on  the  grain  diameter  of  the  sediment  being  considered  and 
on  the  bed  geometry;  therefore,  no  attempt  was  made  to  express  s^  at 
the  bedload  elevation. 

IV.  THi;  SUSPHNDIiD  LOAD  I.\  OSCILLATING  FLOW 
1 . Suspended- Load  llteory  in  Unidirectional  F 1 ow . 


Tlie  suspended- load  theory  for  uiii  directional  flow  and  tlie  available 
field  data  to  test  the  theory  supply  valuable  insight  to  some  problems 
which  exist  in  determining  the  suspended  load  in  oscillating  flow.  For 
this  reason,  Finstein's  (1950)  suspended- load  theor>'  and  field  data  from 
the  Missouri  and  Atchafalaya  Rivers  arc  presented. 

Tl\e  suspension  theory  in  unidirectional  flow  is  based  on  an  equilib- 
rium equation  for  mass  flux  across  a unit  horizontal  area  in  tl\e  flow. 
Assume  the  unit  horizontal  area  is  at  elevation  Y.  Across  this  area 
fluid  is  being  c.xclianged  by  the  vertical  component  of  the  random  motion 
of  fluid  particles  caused  by  turbulence.  From  continuity,  the  picture  of 
fluid  exchange  can  be  simplified  by  assuming  that  tiirough  one-lialf  of  the 
unit  area,  fluid  is  moving  upward  with  an  average  velocity  of  v;  through 
the  ottier  half  area  the  fluid  is  moving  down  with  an  average  velocity  -v. 
If  the  exchange  occurs  over  an  average  distance  of  Ig  it  can  be  assumed 
that  the  downivard-moving  fluid  originates,  as  an  average,  from  an  eleva- 
tion Y + 1/2  Ig  while  the  upward-moving  fluid  originates  from  Y - 1/2  l£3. 
Tlie  important  assumption  is  made  tliat  the  fluid  preserves,  during  its 
exchange,  the  properties  of  the  fluid  at  its  point  of  origin.  If  the 
concentration  of  sediment  at  elevation  Y is  C and  the  sediment  lias 
a settling  velocity  of  Vg , the  equilibrium  equation  for  sediment  flux 
is  given  by; 


C - 4lg  (dC/dY) 
[c  + ilg  (dC/dY) 


j)  (v-VL,) 


- (-v-V..)  = 0 


lliis  ei|uat  i on  reduces  to: 


(-'41 


C Vg  + V (dC/dY)  = 0 . (25) 

To  solve  this  eiiuation  the  term,  1/2  1,  v,  must  be  evaluated.  lliis  is 
normally  done  b>’  equating  this  term  to  the  corros|iond  i ng  term  in  ;i  similar 
equation  of  momentum  exchange;  i.c.,  the  sediment  exchange  coefficient  is 
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assumed  equal  to  the  monentum  exchange  coefficient.  Assuming  that  shear 
due  to  viscosity  may  be  neglected,  compared  with  that  due  to  momentum 
transport,  the  depth,  d,  may  lie  introduced: 

T = [(d-Y)/d]  = -^  p I u - (du/dY)  - [u  + -^Ig  (du/dY)jj  , (26) 

where 

= shear  stress  at  the  bed;  equal  to  gpRS 
T = shear  stress  at  elevation,  Y 
R = hydraulic  radius  (in  feet) 

S = energy  slope 
g = acceleration  due  to  gravity 
p = density  of  the  water 
u = horizontal  flow  velocity 

Using  the  logarithmic  formula  based  on  von  Karman's  (19.34)  similarity  law 
for  the  distribution  of  flovi'  velocity,  du/dY  may  be  calculated: 

du/dY  = (1.0/0. 4)  (u,^/Y)  , (27) 

where  u„  is  the  shear  velocity  and  equal  to  (Tg/p)’^.  Substituting  tliis 
value  into  equation  (26)  and  solving  for  1/2  Ig  v yield: 

jlg  V = (-0.4)  Y u*  (d-Y)/d  . (28) 

Using  this  value  in  e(|uation  (25),  separating  variables  and  introducing 
the  abbreviation: 

Z = V^/(0.4  uj  , (29) 

the  result  can  tie  integrated  from  a to  Y.  The  solution  is: 

(C/C^)  = [(d-Y)  a/Y  (d-a)l^  . (.30) 

It  has  been  found  that  equation  (30)  gives  the  correct  form  of  the 
distribution  function,  but  tlic  value  of  the  exponent  Z given  by  equa- 
tion (29)  docs  not  always  agree  with  the  exponent  tliat  fits  the  measured 
data.  I.et  Z'  be  the  exponent  which  best  fits  tlie  data.  It  was  found 
tliat  for  high  values  of  Z (Z-1.0),  Z'  was  significantly  less.  As  Z is 
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roduci-'d,  tilO  difference  between  Z and  Z'  decreased  and  finall)’  wlien 
Z assumed  values  less  than  unit)',  whicli  is  normally  the  case,  tiie  dif- 
ference between  Z and  Z'  was  small  enough  sucli  that  use  of  equation 
fZyi  allows  accurate  results,  ilie  relationship  Ijetwcen  Z and  Z'  is 
shown  in  figure  28  (f.instein  and  Chien,  1954). 

2 . Si ni lari t i os  Between  Oscillating  and  Unidirectional  Flow. 

'n\o  similarities  between  oscillating  flow  and  unidirectional  flow 
with  low  shear  velocities  are  pronounced.  It  was  found  that  tlie  concen- 
tration distribution  in  oscillating  flow  (from  equation  9)  could  be  ex- 
pressed !)>■ : 


C/C^  = exp(M  Y)  , (31) 

wiiere  M was  found  in  Section  II  to  be  (from  equation  13): 

M = -V3/I:  , (32) 

I:  is  tlic  sediment  exchange  coef'-'icient . Ilie  coefficient,  .'!,  whicli 
defines  the  r.atc  at  which  the  concentration  decays  with  elevation,  beliaves 
in  a manner  similar  to  Z of  the  unidirectional  flow  theory.  In  oscillat- 
ing flow,  :is  in  unidirectional  flow,  the  value  of  M fitting  the  experi- 
mental results  was  elifferent  tiian  tiie  value  which  would  be  predicted  from 
equation  (32),  assuming  F,  indei)endent  of  V..  Figure  15  and  Table  3 
show  that  for  the  four  concetitrat  i on  distribution  curves  obtainct!  when 
V^.  was  increased  from  0.033  to  O.ObZO  foot  per  second,  the  absolute  value 
of  .M  increased,  as  an  average,  by  the  factor  1.19.  Mien  Vg  was  in- 
creased from  0.035  to  0.0498  foot  ])cr  second,  the  absolute  value  of  M 
increased,  as  an  average,  by  1.13.  If  the  exchange  coefficient,  F.,  of 
ct[uat  ion  (32)  were  a function  of  the  flow  hydraulics  only  and  inde]icndent 
of  tlie  sed i ment -sett  1 i ng  velocity,  then  the  average  increase  in  the  abso- 
lute value  of  M would  have  been  0.0620/0.035  (=  1.8)  and  0.0498/0.033 
1=  1.42),  respectively.  Tlierefore,  it  can  be  concluded  tliat  the  sediment 
exchange  coefficient,  F.,  is  a function  of  both  Vg  and  0^,,  and  that 
M,  for  a constant  flow  velocit)’,  is  not  directly  jiroport  ional  to  Vi-. 

This  conclusion  agrees  with  F.instein  and  Cliien  (1954)  that  in  unidirec- 
tional flow  for  high  values  of  settling  velocity  relative  to  the  flow- 
shear  velocity  tlie  sediment  exchange  coefficient  cannot  lie  accurately 
.'ipproxi mated  by  the  momentum  exchange  coefficient. 

Tlie  difference  between  tlic  sediment  and  momentum  excliange  coefficients 
depends  on  the  relative  magnitudes  of  the  sediment-settling  velocity  and 
the  turbulence  intensity.  If  the  sediment-settling  velocity  is  small 
compared  to  tlie  turbulence  intensity,  the  two  coefficients  arc  approxi- 
matel)'  equal;  when  = 0,  the  two  coefficients  are  identical.  As  Vg 
becomes  relatively  larger,  the  difference  between  the  two  coefficients 
incre.ises.  In  unidirectional  flow  the  turlnilence  intensit)’  is  usually 
very  large  compared  to  \(,  ;ind  the  tlieori'  is  accurate  under  most  situa- 
tions. Unfortunately,  this  is  usually  not  tlie  c.ase  in  oscillating  flow. 
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CALCULATED  Z = V*/KU 


Examining  results  b,  c,  and  d of  Section  III,  the  range  of  the  root-mean- 
square  value  of  v'  (s)  is  0.0072  < s < 0.103  foot  per  second.  The 
settling  velocities  of  the  sediments  used  in  the  experiments  were 
Vg  = 0.035,  0.0498,  and  0.0626  foot  per  second.  Vg  is  then,  the  same 
order  of  magnitude  as  the  velocity  scale.  .-\s  discussed  earlier,  the 
sediment  used  in  the  experiments  had  a specific  gravity  of  1.25.  There- 
fore, the  settling  velocity  of  natural  sediment  would  be  even  larger 
compared  to  turbulence  intensity  in  the  approximate  prototype  flow  con- 
ditions of  the  experiments. 

For  illustration,  a typical  oscillatory  flow  condition  of  this  in- 
vestigation will  be  approximated  as  a qu.'is  i -steady  unidirectional  flow 
and  compared  to  the  field  data  of  Figure  28.  Ihis  is  possible  because 
the  time  scale  for  oscillation  is  much  greeter  than  the  time  scale  of 
the  turbulence.  For  example,  the  average  flow  had  a period  of  6 seconds 
and  an  amplitude  of  about  1 foot,  = 0.667  foot  per  second.  A time 
scale  of  turbulence  can  be  defined  as  ^/^o’  5 is  the  thick- 

ness of  the  boundary  layer  and  s^p  is  the  base  vertical  velocity  fluc- 
tuation. If  6 is  defined  as  the  distance  above  the  bed  at  which  the 
boundary  layer  oscillation  velocity  is  equal  to  99  percent  of  the  free- 
streara  velocity  (velocity  given  by  linear  wave  theory  for  y = -d) , 

5 can  be  calculated  from  equation  (5).  This  calculation  indicates  that 

6 is  equal  to  or  less  than  0.05  foot  and  from  Figure  27,  Sq  is  0.11  foot 
per  second.  'ITierefore,  the  time  scale  for  the  turbulence  is  0.45  second 
as  compared  to  a 6-second  time  scale  for  the  oscillation.  To  calculate 
the  theoretical  Z value  for  the  quasi-steady  unidirectional  flow  it  is 
necessary  to  determine  a mean  flow  shear  velocity.  'Hie  flow  shear  veloc- 
ity is  given  by; 

u*  = (g  R S}'^  , (33) 

where  R is  the  hydraulic  radius  (in  feet),  g is  the  acceleration  due 
to  gravity  (in  feet  per  second  squared),  and  .S  is  the  energy  slope  of 
the  quasi-steady  flow.  The  energy  slope  is  obtained  from  Manning's 
equation  by  using  the  root-mean-square  flow  velocity  (=  0.707  2 II  1,/T)  , 
and  estimating  values  of  the  roughness  coefficient,  n,  and  the  hydrau- 
lic radius,  R.  llie  expression  for  the  energy  slope  is: 

s'^  = (n/R°-®^^)  (3  I./T)  . (34) 

Substituting  equation  (34)  into  e(|uation  (33)  yields: 

u*  = (17  n l.)/(T  rO-  167^  _ (35^ 

Using  the  aliove  expression  for  u*  in  ecpiation  (29),  the  2 value  for 
the  quasi-steady  flow  becomes: 


2 = (V,,  T R0-26  7y(,,_g  „ I ■) 


(.36) 
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For  the  average  flow  conditions,  \'g  is  0.035  foot  per  second,  T is 
6 seconds,  and  L is  1.0  foot.  Based  on  the  bottom  roughness  shape  and 
the  liydraulic  radius,  n is  approximately  0.015  foot^’^^^.  'Ilie  hydrau- 
lic radias  is  estimated  from  the  flow  geometry  as  about  equal  to  unity. 
Since  Z is  proportional  to  R to  tlie  0.167  power,  and  therefore  tends 
to  unity,  there  is  probably  not  much  error  in  using  this  estimate.  Using 
these  values  in  equation  (361  yields  Z = 2.1.  In  Figure  28,  this  value 
of  Z is  well  into  the  range  where  the  sediment  exchange  coefficient  is 
significantly  different  from  the  momentum  exchange  coefficient. 

Because  not  enough  data  were  obtained  to  define  a relationship  be- 
tween M and  Vg  and  because  the  sediment  exchange  coefficient  could 
not  be  expressed  as  a function  of  the  shear-stress  distribution,  no 
attempt  was  made  to  derive  a theoretical  relationship  for  the  concentra- 
tion distribution  as  was  done  in  Einstein  and  Qiien  (1954)  for  unidirec- 
tional flow. 

3 . Sediment  Suspens i on  in  an  Oscillating  Flow. 

'Ibis  invest  ig.at  ion  was  done  to  determine  the  behavior  of  sediment 
suspension  in  .a:,  oscillating  flow  and  present  a method  by  which  the  sus- 
pended load  could  be  approximated  from  flow  hydraulics.  It  is  apparent 
from  the  field  measurements  of  unidirectional  flow  (Fig.  28),  and  the 
results  of  this  investigation  for  oscillating  flow  that  the  mechanism  b>' 
which  sediment  is  held  in  suspension  is  complex  and  not  fully  understood. 
For  this  reason,  the  following  method  for  estimating  the  suspended  load 
in  oscillating  flow  as  a function  of  the  flow  hydraulics  is  based  on  the 
general  turbulent  mixing  length  theory  first  proposed  by  O'Brien  (1933). 
His  derivation  is  as  follows: 

Tliere  is  a continuous  up  .•itul  down  motion  of  fluid  across  .any  horiton- 
tal  plane  caused  by  the  turbulent  vertical  velocity  fluctuations.  Ibis 
exchange  motion  is  capal'le  of  transporting  saspended  matter.  Consider  a 
horizontal  reference  section  of  unit  area  :it  a distance  Y from  the  bed. 
The  transfer  of  sediment  in  the  vertical  direction  from  the  region  of 
high  concentration  to  a region  of  low  concentration  through  tliis  unit 
section  will  be  -1/2  1^  v dC/ JY,  where  I^,  is  the  mixing  length  for  the 
sediment  exchange,  v denotes  the  exchange  (.lischarge  through  tlie  unit 
area  due  to  vertical  velocity  fluctuations,  ,iiul  C is  the  concent  r, at  ion 
of  suspended  sediment  witli  settling  velocity,  VV,  , at  elevation  Y.  How- 
ever, a continuous  settling  of  particles  througli  the  unit  are.a  at  a rate 
of  C Vf  exists.  A statistical  equilibrium  condition  is  given  by  equa- 
tion (25).  Tliis  equation  is  identic.al  to  tlie  eipiation  derived  1)>’  Einstein 
(1950),  but  without  the  assumption  tliat  the  origin  of  tlie  sediment  is  the 
s;ime  as  the  origin  of  the  fluid  .and  witliout  assuming  any  dist  ril'ut  ion  of 
fluid  exchange. 

Hie  mixing  length  theory  incorporates  all  of  the  factors  whicli  .affect 
sediuent  suspension  in  two  artificial  v.ariables,  the  velocity  sc.alo  and 
tlie  length  scale.  Ilie  velocity  sc.ale  for  osci  1 l.at  iti);  flow  was  measured 
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directly  in  tliis  investigation  and  is  given  by  equation  (22).  Ttie  sedi- 
ment exchange  coefficient  which  is  the  product  of  the  length  and  velocity 
scales  was  also  measured  directly  and  is  given  by  equation  (13).  From 
those  two  equations,  tlie  corresponding  length  scale  can  be  calculated  as: 

ylg  = (-Vg/M  s^)  exp(-A  Y)  (37) 

For  the  limited  flow  conditions  investigated,  all  the  variables  of 
equations  (13),  (22),  and  (37)  have  been  found  as  a function  of  , 
where  is  a function  of  the  surface  wavelength,  period,  and  water 

depth.  The  variable  A was  found  to  be  a constant  (=  -10.57  feef^), 
and  M and  s^  arc  graphically  given  in  Figures  6 and  27,  respectively. 
Iliese  expressions  for  the  velocity  and  length  scales  are  physically 
reasonable.  .N'ot  only  is  the  base  turbulent  velocity  a function  of  the 
surface  wave  intensity  (Fig.  27),  but  the  exchange  length  tends  to  small 
values  as  the  ocean  bed  is  approached,  indicating  that  no  sediment  should 
be  exchanged  across  the  lied  surface. 

4.  Tlie  Base  Concentration,  . 


With  the  exception  of  the  base  concentration,  , all  the  variables 
needed  to  describe  the  suspended  load  as  a function  of  flow  hydraulics 
have  been  discussed.  'Ilie  following  is  only  a brief  discussion  of  the 
base  concentration.  Kalkanis  (1964)  provides  a complete  mathematical 
derivation. 

In  oscillating  flow,  as  in  unidirectional  flow,  sediment  transport 
is  by  two  different  types;  (a)  Bedload  transport,  and  (b)  suspended- load 
transport.  As  discussed  in  Section  I,  the  thickness  of  the  bedload  layer 
is  about  two-grain  diameters.  Thcrefoi'c,  for  botli  prototype  and  experi- 
mental sediments,  the  bedload  is  contained  in  the  boundary  layer  described 
in  Section  II.  'nie  theory  proposed  by  Kalkanis  to  predict  the  amount  of 
l)edluail  transport  and  tlie  concentration  of  sediment  is  general  and  only 
requires  knowledge  of  tlie  surface  wave  characteristics,  the  water  depth, 
the  bed  seiliment  characteristics,  and  statistical  parameters  which  have 
been  found  experimentally.  Because  the  thickness  of  the  bedload  layer 
is  small,  the  concentration  of  sediment  in  this  layer  is  assumed  const.uit 
and  ei|ual  to  Ibe  concentration,  Co,  for  the  bedload  is  the  base 

concentration  to  be  used  for  the  sus]ieiuled  load. 

Using  as  the  suspended- load  concentration  incorporates  a small 

error  in  the  total  suspended  loail.  As  indicated  previously,  the  bedload 
is  contained  in  the  boundary  layer.  Tlie  distribution  of  sediment  concen- 
tration in  this  area  is  unknown.  Because  the  distance  between  the  top  of 
the  boundar)’  layer  and  the  top  of  the  bed  layer  is  small,  extension  of 
the  exponential  suspension  distribution  ilown  to  the  bed  l.ayer  would  incur 
only  a minor  error  in  the  total  amount  of  seiliment  in  motion. 
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5 . Net  Transport  of  Sediment  in  the  Ocean. 

The  initiation  of  sediment  movement  in  the  ocean  is  by  wave  action. 
Since  the  waves  are  approximately  linear,  the  wave-induced  fluid  motion 
is  a symmetric  oscillation  causing  an  equally  symmetric  movement  of  the 
sediment.  This  motion  ordinarily  cannot  cause  net  transport  of  sediment, 
but  it  does  suspend  sediment  so  that  currents  superimposed  on  the  oscil- 
lating velocity  will  cause  a net  transport  of  the  sediment.  Kxamples  of 
unidirectional  currents  in  the  ocean  are  the  longshore  current  caused  by 
waves  attacking  the  coastline  at  an  acute  angle  and  the  secondary  currents 
caused  by  the  coastline  geometry.  These  two  examples  indicate  that  deter- 
mination of  the  transporting  current  should  be  by  field  measure'.ent . 

6.  Additional  Investigations  .Needed  to  Complete  the  Suspension  Hieor)’. 

The  suspended- load  theory  is  by  no  means  complete.  Tlie  relationship 
between  the  rate  of  decay  of  sediment  concentration,  M,  and  sediment- 
settling  velocity,  \'s  , is  needed.  Figure  15  indicates  that  settling 
velocity  has  a significant  affect  on  the  distribution  of  sediment  con- 
centration, but  there  are  not  enougli  data  to  indicate  the  relationship. 
Additional  measurements  are  required  to  quantitatively  determine  the 
relationship  between  M and  Vg  for  a constant  U^. 

Another  variable  affecting  the  suspended- 1 oad  theory  and  not  studied 
in  this  investigation  is  the  bed  roughness.  It  is  probable  that  a change 
in  bed  roughness  would  affect  the  intensity  of  turbulence  at  a given 
elevation.  Also,  as  discussed  in  Section  II,  tlie  amplitude  of  oscilla- 
tion relative  to  the  wavelength  of  the  bed  roughness  lias  an  affect  on  the 
distribution  of  suspended  sediment.  determination  of  how  the  lied  rough- 
ness affects  sediment  suspension  requires  a great  de.al  of  experimentation; 
however,  the  results  of  this  investigation  will  supply  some  guidelines 
which  would  reduce  the  amount  of  experimental  work  required. 

7 . Cone  1 us  ions . 

These  experiments  in  an  oscillating  flow  simulating  wave  motion  at 
the  ocean  floor  provide  results  which,  in  some  cases,  substantiate  pre- 
vious results  and  provide  new  information  on  the  behavior  of  sediment 
suspension.  It  should  be  stressed  that  the  conclusions  of  this  investi- 
gation are  confined  to  the  rather  limited  range  of  the  varialiles  studied. 

The  conclusions  obtained  from  sediment  concentration  measurements  are; 

a.  'ITie  relationship  between  the  mean  sediment  concentration  and 
elevation  above  the  bed  is  exponential.  Ibis  conclusion  is  based  on  (i.i 
concentration  distribution  relationships  covering  a wide  range  of  proto- 
type flow  conditions  and  each  composed  of  numerous  point  concentration 
measurements.  Typical  examples  of  this  relationship  are  shown  in  Figure 
5 and  the  basic  data  substantiating  tliis  conclusion  are  given  in  the 
appendix.  Although  this  relationship  was  determined  using  an  artificial 
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sediment  i.'ith  a settling  velocity  less  than  found  in  tlie  prototype, 
results  of  experiments  by  otheT  investigators  using  prototype  sediments 
resulted  in  the  sane  exponential  relationship.  In  addition,  the  bed 
roughness  used  by  other  investigators  varied;  therefore,  this  conclusion 
does  not  appear  to  be  limited  to  the  single-hod  roughness  used  in  this 
investigat ion. 

b.  Using  the  above  conclusion  and  the  O'Brien's  (1933)  equation  for 
continuity  of  sediment  exchange  results  in  a sediment  exchange  coefficient 
which  is  independent  of  elevation  above  the  bed.  Behavior  of  the  rate  of 
sediment  concentration  decay  with  elevation,  which  is  related  to  the  sedi- 
ment exchange  coefficient  by  equation  (13),  was  found  to  be  a function  of 
the  flow  velocity  causing  the  suspension  and  the  settling  velocity  of  the 
sediment.  A linear  relationship  between  the  flow  velocity  and  the  sedi- 
ment concentration  decay  rate  was  found  for  the  constant  sediment-settling 
velocity  used  in  the  majority  of  the  experiments.  Tliis  relationship  is 
sliown  in  Figure  6.  Tlie  relationship  between  flow  velocity  and  the  sedi- 
ment concentration  decay  rate  for  other  sediments  studied  in  this  investi- 
gation is  shown  in  Figure  15.  The  limited  data  indicate  a possible  linear 
relationship  for  the  different  sediment-settling  velocities.  Tliere  is  not 
enough  data  to  determine  the  relationship  between  the  concentration  decay 
rate  and  settling  velocity  for  a constant  flow  velocity.  Only  qualita- 
tive conclusions  can  be  obtained  from  the  eight  concentration  distribution 
measurements  shown  in  Figure  15.  For  a constant  bed  roughness  and  flow- 
velocity,  a liigher  sediment-settling  velocity  results  in  a higher  sedi- 
ment concentration  decay  rate.  Ihe  limited  data  consistently  indicate 
that  the  concentration  decay  rate  is  not  proportional  to  the  settling 
velocity  to  the  first  power;  i.e.,  they  are  not  directly  proportional. 

Tliis  and  equation  (13)  imply  that  the  settling  velocity  is  an  important 
variable  influencing  tlie  sediment  exchange  coefficient.  Tlierefore,  in 
oscillating  flows  the  sediment  exchange  coefficient  cannot  be  accurately 
approximated  by  the  momentum  excliange  coefficient  as  is  commonly  done 
in  unidirectional  flow  analysis.  !Jo  experiments  were  conducted  in  this 
investigation  to  determine  how  tiie  afiove  relationships  would  change  witli 
a change  in  bed  roughness. 

Tl\e  conclusions  obtained  from  measurements  of  the  turbulent  velocity 
fluctuations  are: 

a.  The  distribution  of  turl)ulent  velocity  fluctuations  at  a constant 
elevation  in  an  oscillating  flow  was  found  to  be  approximately  normal  with 
a mean  of  zero.  This  relationship  was  determined  from  distril)ut ion  analy- 
ses of  measurements  made  at  two  elevations  above  tlie  bed,  botli  of  which 
were  above  the  boundary  layer  described  in  Section  II.  Results  of  these 
analyses  are  shown  in  Figure  22. 

b.  llie  relationsliip  between  tlie  root-mean-square  turbulent  velocity 
fluctuation  and  elevation  above  the  bed  was  found  to  be  exponential.  Tliis 
conclusion  is  based  on  measurements  of  distributions  made  for  four  differ- 
ent flow  velocities,  all  using  an  amplitude  of  oscillation  of  0.925  foot 
and  covering  apjiroximately  the  s.ime  range  of  prototype  flow  velocities  as 
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studied  in  the  concentration  measurements.  Results  of  these  measurements 
are  shown  in  Figures  23,  24,  25,  and  26,  and  are  tabulated  in  the  appen- 
dix, Table  A-5.  This  relationship  was  valid  for  elevations  of  approxi- 
mately 0.04  foot  above  the  crest  of  the  bed  dunes.  Tlie  relationship  below 
this  elevation,  which  would  be  in  the  boundary  layer,  was  not  determined. 
The  turbulent  velocity  fluctuation  distribution  was  only  measured  for  the 
single-bed  roughness  described  in  Section  II. 

c.  Based  on  the  four  turbulent  velocity  fluctuation  distributions 
described  above,  it  was  concluded  that  the  rate  of  turbulent  velocity 
decay  with  elevation  above  the  bed  is  independent  of  both  the  elevation 
and  the  flow  velocity  generating  the  turbulence.  The  exponential  decay 
rate,  determined  from  a least  squares  curve  fitting  of  the  data,  for  the 
four  distributions  ranged  from  -10.38  to  -10.86  feet'^,  with  a mean  of 
-10.57  feet' and  a variance  of  0.05  foot'^. 

d.  The  relationship  between  the  flow  velocity  and  the  root-mean- 
square  turbulent  velocity  fluctuation  at  zero  elevation  (calculated  from 
the  empirical  relationships)  is  shown  in  Figure  27.  A linear  relation- 
ship is  indicated.  However,  this  relationship,  which  is  far  from  con- 
clusive, is  based  on  only  four  data  points  with  a significant  amount  of 
scatter.  The  qualitative  conclusion  that  the  turbulence  intensity  at 
zero  elevation  becomes  larger  with  greater  flow  velocities  is  not  only 
indicated  by  the  data  but  is  logical. 
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Table  A-2.  Basic  data,  concentration  distribution  measurements,  • 0.035  foot  per  second,  amplitude  < 0.693  foot 


Basic  data,  concdntration  distribution  neasurements , Vg  > 0.0498  foot  per  second,  amplitude  • 0.848  foot 


